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Abstract. We consider particles in M'^, cf > 2 interacting via attractive pair and repulsive four-body 

potentials of the Kac type. Perturbing about mean field theory, valid when the interaction range becomes 

- - - infinite, we prove rigorously the existence of a liquid-gas phase transition, when the interaction range is 

OO ' 

Q.^ ' finite but long compared to the interparticle spacing for a range of temperature. 
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1. Introduction 

\r\ ' An outstanding problem in equilibrium statistical mechanics is to derive rigorously 

^ ■ the existence of a liquid-vapor phase transition in a continuum particle system. While such 
Q"^ ! transitions are observed in all types of macroscopic systems there is at present no proof 
O \ from first principles of their existence in particles interacting with any kind of reasonable 
\ potential, say Lennard- Jones or hard core plus attractive square well. Pair potentials of 
this type obtained by comparison of experiment with low density expansions, accurately 
describe the behavior of both gases and liquids in the ranges of temperatures and pressures 
where boiling and condensation takes place. In fact computer simulations using classical 



\ statistical mechanics of such systems, containing several hundred to several hundred thou- 
Q ' sand particles, convincingly show, via extrapolations which take into account finite size 
O ■ effects, that systems with these type of interactions will have true liquid- vapor phase tran- 
sitions in the thermodynamic (infinite volume) limit. In this paper we go some way toward 
^ ! a proof of such a transition, i.e. we prove for the first time the existence of a liquid vapor 
^ \ transition in a continuum particle system with finite range interactions and no special 
symmetry. 

Historically the first proof of liquid- vapor type phase transitions was given for lattice 
systems (which are isomorphic to Ising spins) . These systems can be thought of as idealiza- 
tions of real fiuids in which however the natural continuous spatial translation invariance 
symmetry is replaced by that of the lattice Z^, (i > 2. It was Peierls [P] who first gave 
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a convincing argument (later made fully rigorous by Dobrushin [Dl] and Griffiths [Gr]) 
of the coexistence of different phases in such systems. The power of Gibbsian statistical 
mechanics to produce such rigorous results was brought home to the general science com- 
munity by the dramatic work of Onsager [O] explicitly solving the two dimensional Ising 
model (or lattice gas), with nearest neighbor interactions. Since that time there have been 
found many other exactly solvable two dimensional lattice systems [Bax]. At the same 
time the development of various types of inequalities as well as the powerful Pirogov-Sinai 
formalism [PS] have resulted in a comprehensive rigorous theory of phase transitions in 
lattice systems, at sufficiently low temperatures. 

For continuum systems there is at present proof of a phase transition for the two 
component Widom-Rowlinson model [WR] in which the dominant interaction is a strong 
repulsion between particles of different species. The physical arguments of WR were made 
rigorous by Ruelle [Rl] who was able to generalize the Peierls argument to prove phase 
coexistence in this system. Ruelle's proof strongly exploits the symmetry between the two 
components present in the WR model. The same is true, at least to some extent, of various 
extensions of this model [LL], [BKL], [Ge]. For continuum systems without some special 
symmetry the only proofs of phase transitions so far are for systems with interactions 
which decay very slowly or not at all. Such one dimensional models with many particle 
interactions were analyzed and proven to exhibit phase transitions by Fisher and Felderhof 
[FF] . More recently Johansson [J] has considered interactions in one dimension which decay 
as r~", q; e (1, 2), proving that at low temperatures there is phase transition in the sense 
that the pressure is not differentiable. 

The mean field or van der Waals type of phase transition was also first derived rigor- 
ously by Kac, Uhlenbeck and Hemmer [KUH] for a class of one dimensional models, hard 
rods of radius one with an added pair potential 

(l>'yiQi,Qj) = -Oi^l^w[-l\Qi- Qj\], l,a>0 (1.1) 

in the limit 7 — » 0, see also Baker, [Bak]. This was later generalized by Lebowitz and 
Penrose [LP] to d-dimensional systems with suitable short range interactions and Kac 
potentials of the form 

(f>fiQii(lj) = -Q!7'^J(7ki - Qjl) (1-2) 

with 

[ J{r)dr = 1, J(r) > (1.3) 

In the thermodynamic limit followed by the limit 7^0 the Helmholtz free energy a takes 
the form, for a fixed temperature 

lim a{Q, 7) = CEiaoiQ) - \aQ^] (1.4) 

Here q is the particle density, ao is the free energy density of the reference system, i.e. 
the system with a = in (1.2). ao is convex in q (by general theorems) and CE{f{x)} is 
the largest convex lower bound of /. For a large enough the term in the curly brackets 



in (1.4) has a double well shape and the CE corresponds to the Gibbs double tangent 
construction. This is equivalent to Maxwell's equal area rule applied to a van der Waals' 
type equation of state where it gives the coexistence of liquid and vapor phases [LP] ; see 
also [vK]. 

In this paper we prove the coexistence of liquid and gas phases for systems with finite 
range interactions as small perturbations at finite 7 > from the mean field behavior at 
7 = 0. The same approach was recently applied to the lattice case for Ising models, [CP], 
[BZ], [BP], where the Peierls argument applies directly, because of the spin flip symmetry 
of the model. The absence of symmetries in our case requires instead the whole machinery 
of the Pirogov-Sinai theory. This theory was developed for and is very powerful in proving 
phase transitions of lattice systems at low temperatures, when the pure phases are close to 
the ground states. To apply the theory to our system we will perform a coarse graining in 
which we divide the space into large cubes of side £ and introduce variables Qx which are the 
particle densities in the cubes labeled by x. After integrating out all the other variables, 
we will be left with a new system described by the variables Qx- We will then show that 
the distribution of the is still Gibbsian with a new Hamiltonian and temperature. The 
essential point is that the ratio between the new and old temperatures scales as By 
suitably choosing the side of the cubes, we will then enter into the low temperature regime 
where the Peierls and the Pirogov-Sinai (hereafter denoted by P-S) methods apply. In this 
new perturbative scheme, the unperturbed state is described by mean field (formally 7 = 0) 
and the small parameter of the expansion is the inverse interaction range 7, instead of the 
temperature in the traditional approach. By choosing a suitable range of values of chemical 
potential and temperature we will then be able to put ourselves at the vapor-liquid phase 
coexistence. 

To insure stabilization against collapse, which would be induced by a purely attractive 
pair potential, the natural choice is to replace the point particles by hard spheres. Our 
approach however does not work in such a case, as we need a cluster expansion for the 
unperturbed reference system (i.e. without the Kac interaction) at values of the chemical 
potential or density for which it is not proved to hold. We avoid the problem by considering 
point particles and insuring stability by introducing a positive four body potential of 
the same range as the negative two body potential. In this way we avoid having to 
control strong short range interactions, something beyond our present day abilities for 
dense continuum systems. 

2. Definitions and Results. 

For ease of reference we give below the main definitions followed by the main result. 

Particle Configurations and Phase Space 

We consider a system of identical particles in M*^, d > 2. Particles are points in M'', 
particle configurations are countable, locally finite collections of particles. The set of all 
particle configurations in A, A C M.^, is the phase space Q^^^ and we simply write Q when 
A = W''. The particle configurations are denoted by q or by q^^"* when we want to underline 
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that they are in Q^^\ We write = (gi, .., g^) when the configuration consists oi \q\ = n 
particles positioned at points gi, G 

Free Measure 

For a bounded measurable subset A of the free measure dq on Q^^^ (also called 
the Liouville measure) is 

/ dQ f{Q) = y2-i I dqi--- I dqn f{qi, ..,qn) (2.1) 

where / is any bounded measurable function on Q,^^\ 

Hamiltonian 

The energy of the configuration q = (qi) is given for our system by the formal Hamil- 
tonian 

H^Ai) = -Mil -^^Yl 

^ ' " (2.2) 

+ ^1X1X1 Yl Yl A'^^ ^^^^ ' > (lis > (Ha) 

Here A is the chemical potential while ■) and Jj^\-, -r) are respectively two and 

four-body potentials. For notational simplicity we choose 76 {2~", n G N}; 7 is a scaling 
factor in the definition of J^^\-, •) and J^'^\-, •, •, •) which are chosen of the form 

4^'fe.,fe)=7-/*nv,. ,<-,-.«. (2.3) 

r ^ 

J^''Hqn,qi,,qis,Qu) = l''' dr l[lt _ (2.4) 

where I is the characteristic function of the set A and Rd is the radius of the ball in M.^ 

A 

having a unit volume. Denoting by B^{r) a ball of radius ^~^Rd centered at r e IR*^ we 
have 

j\^\qn,qi,) = l^''\B^{qi,)nB^{qi,)\ (2.5) 

and 

4''\qn,(li„qis,qu) = l^'^l n|=i S^(gi,)l (2-6) 

Let J^^^ = j}^^ and J^^^ = Then the scaling properties of j!y^^ and </y'^^ can be 

expressed by 

4^Kqn,qi,) = l^J^^Kiqi^iqi,) (2.7) 



and 

2-5.4 4^^fei,9i2,9i3,9i4) = 7^'^<^^^^(79ii,---,79i4)> (2-8) 

which is similar to (1.2). It is also clear that 



2.5.5 



J J^^\0,r) dr^ J J^^^(0,ri,r2,r3) dridr2dr3 = 1 (2.9) 



as in (1.3). 



As both and are positive we have in (2.2) a competition between an attractive 
pair and a repulsive four-body potential. When the scaling parameter 7 is small (but 
finite), the model has a large but finite interaction radius, 2Rdj~^, and a small interaction 
strength between any given two or four particles. Nevertheless, the total strength of the 
interaction between a given particle and all other particles in a configuration q of bounded 
nonvanishing density is of order 1. These are characteristic properties of the Kac potentials 
which, as noted earlier, usually reproduce the van der Waals theory [LP] in the scaling limit 
7^0. The specific form (2.3)-(2.4) of the interaction J^'^^ and J*^'^-' makes the analysis 
simpler, see the discussion in the beginning of Section 3. 

Gihhs Measures 

For any two configurations q = (qi) and q = (qj) denote by g U g = {qi,qj) the 
configuration containing all particles from both q and q. The conditional energy of q^^^ 
given a configuration g e Q is 

2.3 if^,A(g(^)|g) = H^Aq^""^ U g) - H^^Q) (2-10) 

This conditional energy consists of two parts: the energy, Hj^x{q^^^), of the configuration 
g(^^ itself and the interaction energy 

2.3.1 [/^(g(^)|g) = H^Aq^^^Iq) - H-yAi^""^) (2-11) 

between configurations g*-^-* and g. 

Let A"^ = \ A be the complement of A. Given an inverse temperature (3 the Gihhs 

measure fi^^^ ^{dq^^^q^^ )) in a bounded measurable set A with a boundary condition g(^ ) 

is the following probability measure on Q^^^ 



where S..y^^^;^(A|g(^'^)) is the partition function 



p A(A|g-(^^)) = / dg(^)e-^^--(^^^^l^"'^ (2.13) 

The infinite volume Gihhs measures fXj^i3^\{dq) are probabilities on Q such that for 
any bounded measurable set A and /U^,/3,A-alniost any boundary condition g*^^ ^ G Q^^ ■* 

the conditional measure A*7,/3,a(<^'?*-^'*|^*-^ ■*) is equal to /U^^^ ^^{dq''^'^ \q^^ )) given by (2.12). 

We say that a translation invariant Gibbs measure has a particle density ^ > if for 
any bounded set A the expectation of |g fl A| is equal to g\A\. 
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The Main Result 

Our main purpose is to investigate the phase diagram of the model (2.2) and to 
prove that a phase transition of the hquid - vapor type takes place for some values of the 

3 

parameters. Let Pc = (f ) ^ and Pq > /?c be a number defined via (3.13) below. 

s2.i Theorem 2.1 

For any P G (/3c, Po) there exist functions 'foiP) ^nd A(7, P) such that for < 7 < 7o(/?) 
the model (2.2) has at least two distinct Gibbs measures ^ These measures 

are translation invariant and ergodic, with an exponential decay of correlations. They 
have particle densities respectively equal to g* _ > and g* p _^ > _. The quantities 
A(7, P) and p ± have limits as 7 — > 0, denoted by X{P) and gf3,±, and there exist positive 
constants c and a such that, |A(7,/?) — X{P) \ + Yla=± \ (3 a ~ Sf3,cr\ ^ ^7*^. 

The limit quantities X{P) and ^/3,± are given exactly by mean field type formulas 
[LP] and our proof of Theorem 2.1 is a perturbation theory constructed around the mean 
field picture carried out by using the P-S theory [PS]. The restriction P G {Pc, Po) rather 
than P > Pc allows for a simplified proof; we hope to present a proof for all P > Pc in a, 
future paper. While technically different and applied to different settings our methods are 
nevertheless conceptually close to those of [DZ]. 

Contents of the Remaining Sections 

In Section 3 we give an outline of the proof formulating a number of statements which 
are proven in subsequent sections. In Section 4 we prove Peierls estimates on contours, and 
in Section 5 we use cluster expansion to investigate properties of the effective Hamiltonian 
obtained after the coarse graining transformation. In Section 6 we study the restricted 
ensembles proving the basic property of the P-S scheme, namely that it is possible to 
adjust the value of the chemical potential in such a way that the pressures in the two 
restricted ensembles are equal. We will conclude from this the proof of Theorem 2.1. Two 
auxiliary statements used in the proof are included in the text as Sections 7 and 8. 

3. Scheme of Proof. 

The starting point of our approach, which as explained in the Introduction is of 
crucial importance, consists of coarse graining: we partition the space into cubes ct'^\ x 
the centers of the cubes, of side small positive number, (the subscript 2 

foresees the use of other scales that will be introduced later in the proof). Given a particle 
configuration q, we call g^;, the particle densities in each cube and we integrate out all 
other variables, i.e. we consider the marginal over the variables g^- The new measure is 
still Gibbsian and its new, effective Hamiltonian (written as a function of the gx) can be 
characterized with a remarkable accuracy by cluster expansion techniques, see Section 5. 
As mentioned earlier, the main point of this procedure is that the new effective inverse 
temperature becomes P£2- We can thus enter into the very low temperature regime by 
taking 7 small enough, with £2 correspondingly large (in Section 5 we will absorb the 
temperature into the Hamiltonian). We are then in the right setup for the P-S theory. 
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2.13 



An analysis a la [LP] , which we omit here, would show that the new effective Hamiltonian 

converges formally in the limit 7 — ^ to the mean field free energy functional J^p^xig); see 
(3.1) below. We will begin our analysis by characterizing the ground states of J-p^x^g) and 
thereafter use the P-S theory to investigate the perturbations at 7 > 0. 

Mean Field Free Energy Functional and Ground states 

The mean field Gibbs free energy functional J-p^xio)-, Q = Q{r) denoting a non- negative 
bounded measurable function in R^, is 

^/3,a(^) = / dT^{\ogQ{T)-l)- I drXg{r) 

dridr2j'''^\ri,r2)g{ri)g{r2) (3.1) 

+ dri---dr4j'^^\ri,...,r4)g{ri)---g{r4) 

The first integral is the entropy contribution to the free energy (more precisely the product 
of the temperature times the entropy of the ideal gas with the changed sign). The other 
three terms arise from the corresponding interaction terms in (2.2). More details concerning 
the relation between (3.1) and (2.2) can be found in Section 4. 

The mean field ground states are, by definition, the minimizers of Tp^xio)- To find 
them we set 

n{r,g)= I dngin) (3.2) 

J\r-ri\<Rd 

and 

T(r, g)= [ dr^^{ log g{r,) - l) (3.3) 

J\r-rr\<Ra P 

TZ{r, g) and X(r, g) are respectively the mean density and the mean entropy of g in the ball 
B{r). With this notation we can rewrite (3.1) as 

2.15 J^p,x{q) = Jdr (x{r, g) - Xn{r, g) - ^n{r, gf + ^7^(^, g)"^^ , (3.4) 

which is true only because of the special form, (2.3) and (2.4), of J^^^ and ,7^^). In fact 
(3.4) is the main reason for choosing J^^) and J^^^ of this form. By convexity 



2.14 



2.14.1 



2.16 



I{r,g)>^^^[lognir,g)-l) (3.5) 



and equality is achieved only if 7^(r, ^) = g{f). 

It follows from (3.4) that if s* > is a minimizer of the function 

2 4 

2.17 F^,x (.) = ( log . - 1) - A. - |j + ^ (3.6) 
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then Q{r) = TZ{r, g) = s* is a minimizer of J^i3,\{q). The second derivative 



2.20.4 



2.22 



2-17.1 Fl^{s) = -1 - 1 + ^ (3.7) 



of Ffj^x{s) is positive ii (3 < (3c — (|) ^ . Hence for (3 < (3c and any A the function Fg^^l-s) 
is convex and has a unique minimizer which is the root of the equation 

2.20 -s + j^ s^ + ^ log s - A = (3.8) 

For (3 > (3c this equation has three positive roots with two of them, s = Qf3,\,- and 
s = Qf^^x^^, being local minimizers of Fg^^l-s)- Furthermore, there exists a unique A = A(/3) 
for which both local minima are the global ones and the function -F/3,a('S) has a "double 
well" shape with the same "depth" of the wells. We set g/3^± — Qp,x{i3),±- Clearly for 
(3 > (3c and A = A(/3) the densities Q{r) = g^^- and g{r) = gp^s^ are distinct mean field 
gound states. For later purposes we remark that 

2.20.1 -1 < \{(3) < (3.9) 
which can be checked by direct calculation. 

Mean Field Equations and Contraction Property 

As the ground states are minimizers of J-'/3,\{g), they satisfy the mean field equation 
^^P,\{q)/^q{^) — 0. By an explicit calculation we then find that they are fixed points of 
the transformation 

2.20.2 Q{r) ^^{g{-),r) (3.10) 
where 

^{g{-),r) :=exp|A+ J driJ^^\r,ri)g{ri)- J dn.. J drsJ'^^\r, ..,rs)g{ri)g{r2)g{r3)^ 

2.20.3 (3.11) 
Setting g* = gp_ or g* = + the derivative ip{r) = S^(^g{-),r)/Sg{r) computed at 
g{r) = g* is equal to 

i;{r) = (3g-J^'\0,r)[l-{g*f/2] (3.12) 



We define Po in Theorem 2.1 so that for all /? e {(3ci(3q) 

J dr\tP{r)\ = pg*\l - g'''^ /2\ < 1 (3.13) 

If (3.13) holds, the transformation $ is a contraction (in a neighborhood of the ground 
states and using sup norms). 



2.20.5 



The existence of Po > (5c follows by noting that since (5q* <1 (by (3.7) and because 
^p,xiQ*) > 0) the only condition to check for (3.13) to hold is Pg*[l - Q*'^/2] > -1. By 
(3.7) it can be rewritten as 

F^J<^*) < ^ (3.14) 

and the result then follows from the fact that ^{g*) — > as /3 — > /3c while both ^^3^-1- 
remain bounded. 

We will see in Section 5 that for (3 G (/3c, Po) and 7 > small enough, the Dobrushin 
uniqueness condition is satisfied by the effective Hamiltonian of the system restricted to the 
ground state ensemble (defined later in this section). This uniqueness condition appears 
to be the 7 > analogue of the contraction property (3.13). We will in the sequel restrict 
to P e (PcPo) which guarantees the following properties. 

For Pc < P < Po there exists a positive number ({P) such that 



2.21 



max P{gi3,at + 2C(/3)) 

cri=±l,cr2 = ±l 



1 - 



a{p) < 1 



(3.15) 



Consequently there exist a positive number 5{P) such that for any A e (A(/3) — 5(/3), A(/3) + 
5{P)) one has 



2.21.1 



max P{gi3,\,a^+C{P)) 

CTi=±l,a2=±l 



<a{p) 



(3.16) 



We expect that for 7 small enough the particle configurations will have densities close 
to the mean field ground states. To investigate the issue we need a preliminary spatial 
scale for computing particle densities and a notion of closeness between densities. We will 
then have a picture of the particle configurations in terms of spatial regions where there 
is agreement or disagreement with the ground states. 

The Partitions T>^^\ the Cubes Ci^^ and the Densities gx^ 

Let 'D^^\ £ e {2", n G N}, be decreasing partitions of into cubes C^^^ of side i.e. 

is coarser than if £ > £'. For any r G R'* we denote by C^^\r) the cube of "D*^^^ 
that contains r. We suppose that the centers of the cubes C^^^ are in Z^. Consequently 
the centers of cubes C^^) are in = ^Z*^ + . . . , ^) . For a; G we denote Ci^^ the 
cube of X>(^) centered at x. 

Given a region A, [A]*^^^ is the maximal V^^^ measurable subset of A, i.e. the union of 
all C^^) contained in A. We also identify [A](^) with {x e Zj\ x e [A](^)}. The R'^ volume 
of [A] is denoted by | [A] \ while its Zf volume, i.e. the number of lattice points, is 
denoted by ||[A]W||. 

The density of a configuration q in a box Ci^^ IS 

(3.17) 



3.2.1 



(g) = t 
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Accordingly g^^\q) = {Qx (q)) is called a (D^^) measurable) density configuration corre- 
sponding to a particle configuration q. We denote ^*^^)(A) = {gx\A)) a density configura- 
tion in [A]*^^-* not related to any particle configuration. 

The notation [q]^^^ is used for the particle configuration obtained from q = (qi) by 
shifting all qi to the centers of the corresponding boxes C'^^\qi). For a configuration q^^^ 
we set 

H,^^{e^'\q^^^)) = H,,xi[q^^Y^) (3.18) 

and 

H.A^^Hq^'^W''^) = H,,,{[q(^^f^q(^'^), (3.19) 
where A is a region and q^^""^ is a boundary condition. 

Spatial Scales, rj Functions and Ground States Configurations 
Several scales are of special interest for us 

3.2 £i= £1,^ = 7-1+"!, £2 = £2^^ = £3^£3^^=^-l-«3^ (3.20) 

where the numbers < ctj < 1 are rational, cci = 1/2 and 0:2 + 0:3 <S {2d)~^. 

The scale £2 is our coarse graining scale. Agreement or disagreement with the ground 
state is indicated by r^-functions, rj = r]{q) = {r]x{Q))j ^ £ '^e^' 

r-1, foT\g^^'\q)-Qp,_\<C 
3-2-2 Vx{q) = < +1, for \g^I'\q) - g/3,+ \ < C (^-^l) 

I 0, otherwise 

( — C{l3) being taken from (3.15). Particle configurations q or density configurations Q^^\q) 
are called compatible with rj ii rj = ri{g^^\q)) = rj{q). 

We say that a particle configuration q belongs to the ground state ensemble of the 
± phase or equivalently liquid/vapor phase if r}{q) = ±1. For brevity we will sometimes 
simply say that g is a ± ground state configuration. 

((■) 

For a scale I we define the corresponding standard ground state configurations q_ and 
q^^ as those which have Q^^\q^l^) = Qp^- and Q^^\q^^^) = with all the particles placed 
at the centers of the corresponding boxes C^^\ When £ = £2 we will drop the superscript 
(£). 

The other scales, £3 and £1, are used to construct contours (see below) and to do some 
approximate calculations (see Section 4) respectively. 

Our proof of the phase transition or coexistence of phases is based on the following 
qualitative picture describing typical particle configurations in the two pure phases. In the 
gas phase a typical configuration q coincides in most of with some typical configuration 
in the gas ground state ensemble. Inside this "sea" of the gas ground state there are 
rare "islands" occupied by the liquid ground state. These two types of ground states 
are separated from each other by regions in IR*^ which are called Peierls contours. The 
"excess" free energy of q with respect to the free energy of the gas ground state occupying 
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all of M.^ is concentrated in the vicinity of these contours and is proportional to their 
volume. A similar, inverse, picture describes typical configurations of the liquid phase; 
the liquid ground state now forms a "sea" with rare "islands" of the gas ground state. 
Thus the typical configurations of the gas and liquid phases are distinct with the density 
of particles being an order parameter distinguishing them. A rigorous verification of the 
above picture is fairly straightforward for simple lattice gases; Ising systems, in which the 
gas and liquid ground states correspond respectively to "all sites empty" and "all sites 
occupied" . It requires however considerable work for continuum systems and we start with 
precise definitions. 

Correct and Incorrect Sets, Contours 

A family of cubes from the partition I>(^) is called -k- connected if its closure is a 
connected subset of M.'^; we consider from now on d > 2. We then say that r is a (+)- 
correct point of a configuration q if ri{q) = +1 in [C^^^\r)]^^'^^ and in each of the neighboring 
25(^3) measurable cubes T*c-connected to C^^^\r). Similarly one defines (—)- correct points. 
Finally, if r is not a correct point of q it is an incorrect one. The connected components 
of the incorrect points of q form the supports of the contours. The pair consisting of the 
support of the contour and the restriction of r]{q) to this support is called a contour of q 
and is denoted T{q). Observe that any q uniquely defines its contours Ti{q). 

Axiomatically a contour T = (Supp(r), 77^) is defined as a pair which consists of a 
bounded ★-connected V^^^^ measurable set Supp(r) called the spatial support of F and 
a spin valued function 77'" on Supp(r), with the condition that there exists at least one 
configuration q that gives rise to F. 

To motivate the next definitions and to outline further constructions let us consider 
first an oversimplified problem. 
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A 




A 



Fig.l 

We choose boundary conditions q^^ ^ in the + phase and we impose that in the region 
A of Fig.l the configuration is in the + phase, in D it is in the — phase while in Supp(r) 
it agrees with rj^ . 

The simphfied contour ensemble gives rise to the simplified partition function Z (we 
are using here the notation that will be discontinued in the sequel) 



z = y I dq^""^ I a 



(3.22) 



Here the external sum is taken over all contours F of the type described by Fig.l and 
the goal now is to rewrite Z as an integral over the + ground state ensemble, i.e. as the 
integral over q'^^^ with {ri{q^^^) = 1}. For B and C strips of width 7"-'^ in Supp(r), as in 
Fig. 1, we write 



Z = Y, dq^""^ a 



where the statistical weight W{T\q^^^) is 

w^(r|g(^)) = 



(3.23) 



(3.24) 



jQ(Supp(r)\B) v{<l)=V Jq(d) 



jQ(Supp(r)\B) ri{q) = l J q^o) viq^ ') = 1 
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In (3.23) the configuration q^^^ is not related to r]^ and the contours F appear only through 
their statistical weights as extra variables in the partition functions. This is the contour 
model in the present over-simplified context. The advantage of the contour model is to 
work in a ground state ensemble, the price is the extra variables F. The game, in the 
general case with many contours, is to get good estimates on their statistical weights 
(3.24) in order to control their proliferation. 

Let us examine this issue still in the context of the oversimplified example. We should 
think of the integral in (3.24) over the region Supp(F) \ S as a "surface term" and of the 
integral over D as a "volume term". The second one looks therefore more likely to take 
large values and we consider it first, as it appears in the numerator in (3.24) denoted by 
N. The integral is over all configurations in the region D that are in the ground state 
ensemble of the — phase. The boundary conditions for this partition function are fixed in 
the region C and they are also configurations of the — phase. In fact the — phase extends 
over the whole £3 cubes ★-connected to D, recall the definition of contours. 

The intuition (which is made rigorous in Section 6) is that since C is well inside the 
— phase, then the configurations in a neighborhood of C arc actually not only in the — 
ground state ensemble but they are very close to the standard ground state configuration. 
We thus write (recall that N denotes the numerator in (3.24)) 

3.8el N = T1T2T3 (3.25) 

where Ti, T2 and T3 are defined as follows. 

Splitting 1/2-1/2 the interaction between Supp(F) and D and recalling the intuition 
that the configurations involved in such interactions are very close to the — ground state 
configuration, we write 



3.8e2 



3.8e3 



3.8el.l 



3.8e4 



7Q(Supp(r)\B) vi<l)=V 

Denoting by /a,- the thermodynamic limit of the pressure associated with the parti- 
tion function restricted to the — ground state ensemble we set 

T2 = e^^-l-^l (3.27) 

The third term T3 is implicitly defined so that (3.25) holds. T3 takes into account the 
errors made by replacing in (3.26) the interaction 

it/,(,(^)|,(^)) - C/,(^(^)ki^^) - lu.iqi'^^ie^) (3.28) 

and the error which comes from T2 where the partition function is replaced by its thermo- 
dynamic limit in the sense of (3.27). 

The error term T3 appears to be bounded by 

I logT3| < C7^/^ jsurface of D\ (3.29) 
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where c is a constant. The above inequahty is proven in Section 6 in the general case. 

After an analogous decomposition in the denominator, we see that the dangerous 
volume terms simplify between numerator and denominator if we are able to choose the 
chemical potential in such a way that /a,- = /a,+- Once again, the possibility to solve 
this equation relies on a rather explicit representation of the pressures, /a,±) which one is 
able to derive. 

The ratio of the terms Ti (from numerator and denumerator) is bounded by 

3.8e5 exp(-c£^£3'^C^|Supp(r)|) (3.30) 

This is the famous Peierls estimate which is proven in the general case in Section 4. Its 
proof is relatively simpler than the proof of (3.29) and relies directly on properties of the 
mean free field energy functional. This step is similar to the analogous one in Ising models 
[BP]. 

Combining (3.29) and (3.30) and recalling that dai + da2 is assumed to be much less 
than 1/4 we conclude that the whole statistical weight W{T\q^^^) is bounded as in (3.30). 
Let us now go from the oversimplified example to the general case. 



More About Contours 
Given a (large enough) region V and £ > we call 

3.3 a^^V = {r e : dist(r, V) < £}; S^^'^V = {r e V : dist(r, V) < £} (3.31) 

and we set dV = d^^~'^V and 5V = S^^~'^V. 

Denoting by Ext(r) and by Intm(r) respectively the unbounded and the bounded 
connected components of (Supp(r))^ we observe that ri^ takes the same values (j{T) and 
cr^(r) on each of the sets d^^^^'Ejxt{T) and d^^'-^^lntm(r) correspondingly. The regions 
Ext(r) and Int(r) = Umlnt^(r) are respectively called the exterior and the interior oi F, 
in the previous example they are A and D. 

For a contour F denote 

3.5.0 5m{T) = aint^(F), (3.32) 
(the set C in the example of Fig. 1) 

3.5 5={T) = aExt(r) U I U 5^{V) j , (3.33) 

\m: a^(r)=<r(r) / 

(the set B in Fig. 1, the last term being absent in Fig. 1) 

3.5.1 5^(r)= U 5m{V), (3.34) 

m: a^(r)7^a(r) 

(this is the set C in Fig. 1) and 

3.5.2 5{T) = 5={T) U 5^{T) = (5Supp(F) (3.35) 
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(this is S U C in Fig. 1) 

Two contours are called disjoint if their supports are not ★-connected. We say that 
a contour belongs to a region if its support is a subset of that region. In any collection 
of pairwise disjoint contours there is a uniquely defined subcoUection of external contours, 
i.e. contours which do not belong to the interior of any other contour in the collection. 

A contour T surrounds a point r if r e Supp(r) U Int(r). A contour T3{q) separates 
contours Tiiq) and T2{q) if Supp(ri(g)) G Int(r3(g)) while Supp(r2(g)) G Ext(r3(g)). 

Consider an arbitrary configuration g which differs from a ground state configuration 
of a given phase a only inside a bounded region A. Then all contours of g belong to A and 
for any external contour r*^'^*(g) one has a{T^^^) = a. Note that non-external contours 
r(g) may have o-(r(g)) = —a. Moreover, the contours of g satisfy the following matching 
condition. If Supp(ri(g)) G Intm,(r2(g)) and there is no r3(g) separating ri(g) and r2(g) 
then cr(ri(g)) = a"m,(r2(g)). This matching condition is highly non local and represents 
the main difficulty in understanding contour statistics. 

The idea of the P-S theory is to remove matching condition and to construct an 
equivalent contour model with no matching rules present. Generally this can be done by 
modification of the statistical weights of contours and shifting the difficulty to the estimate 
of these modified statistical weights. 

Generally the study of the distribution of contours is important because to prove the 
existence of distinct translation invariant + and — phases it is enough to show that the 
probability of the event that an external contour surrounds the origin is sufficiently small. 
In a translation invariant measure this implies that only finitely many contours surround 
any given point r. Then one of the ground states occupies the infinite region njExt(r|''*(g)) 
which is exactly the "sea" discussed earlier. Furthermore, the probability of the event that 
?7o(g) = 1 (here denotes the origin) can be taken as the order parameter distinguishing 
two phases. This probability is less than 1/2 for the — phase and it is larger thanl/2 for 
the -|- phase as follows from a bound on the probability of a contour of the form (3.30) 
(which will be shown to be true for all contours). As we will see even in the general case 
when there is more than just one single contour in the contour model the bound (3.30) 
implies that the contours are very rare and one can control them with an analysis not too 
distant from the oversimplified example. 

Let us recall now that the bound (3.30) on the statistical weight of a contour was 
obtained after adjusting the chemical potential so that the pressures in the two ground 
states ensembles are equal. In the general case this is not a simple task because there are 
contours inside contours, i.e. inside the regions of type D of the example in Fig. 1 there 
are other contours and so on. One needs to treat contours recursively as in the one-contour 
example to obtain at the end of the recursion the contour model with the properly modified 
statistical weights of contours. Only after that may one try to equalize pressures in these 
contour models. The problem is that even if the bound (3.30) holds the pressures in the 
contour models will depend on the contours and we are then in a sort of loop: to have the 
good bound (3.30) on the contours we need to make the pressures equal, but to control 
the pressures we need a good bound on the contours. 

As explained by the P-S theory it is possible to survive such an impasse, we will do 
it by following the Zahradnik approach [Z] . We will define a new statistical weight with a 
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cut off which by its definition cannot exceed the value (3.30) (with a suitable constant c). 
In this context it will be easy to iterate the procedure of the example with one contour 
to derive a contour model representation of the true partition function. Such auxiliary 
partition functions (one for each phase) give rise to the corresponding pressures and we 
will adjust the chemical potential so that these two pressures (one for each phase) are 
equal. It will then turn out (see [Z]) that for this particular value of the chemical potential 
the contours satisfy the bound (3.30) without the need of the cutoff so that for this special 
value of the chemical potential the auxiliary and true partition functions are equal. 

The next subsection contains exact definitions and statements which are necessary for 
application of general P-S theorem (see [PS], [Z]). 

Auxiliary Partition Functions 

Given a phase a and boundary condition q^^ ^ belonging to the ground state ensemble 
of this phase we define the auxiliary partition functions and the truncated statistical weight 
as 

and 

where c = C(3.3o) 

^A^P|-(5=(r))) (3.38) 



/ dqTL e-^^-^(^l^"''^'"") ff Z^. . (lnt^(r)|?(^i-t-(r))) 



7e(Supp(r)\.= (r)) r,(«)=<7(r) 11 7,/3,A V -V J 

To show that the definition is well posed, we first say that a statement related to a V^^^ 
measurable region is proved by induction in volume if this statement is true for a V^^^ 
measurable region A as soon as it is true for all V^^^ measurable subsets of A. We then 
observe that (3.36)-(3.38) constitute a single inductive definition. At the initial step of 
this induction in volume one considers only contours with empty interior and uses (3.37)- 
(3.38) to define their truncated statistical weights. Then (3.36) allows one to calculate 
the auxiliary partition functions for sufficiently small regions admitting inside them only 
contours without interior. Afterwards (3.37)-(3.38) are used again to define the truncated 
statistical weights of contours having sufficiently small but non empty interiors and so on. 

We hope that a notational similarity of the numerator and the denominator of (3.38) 
does not hide from the reader the fact that these expressions are very different. For 
example, in the numerator of (3.38) configurations qi9'i-n^m{T)) belong to the ground state 
ensemble of the phase a = o"(r) while in the denominator of (3.38) these configurations 
belong to the ground state ensemble of the opposite phase, —a. 
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s3.2 Statement 3.2 

If the truncated statistical weight W'^{T\q^^ i'^))'^ of any conto^ir is smaller than the 
second term in the argument of the min in (3.37) then the auxiliary partition function 
coincides with the true one. 

The proof of the statement is standard in the P-S theory and we omit it (see [Z]). 

Another standard observation is that for any region A with the boundary condition 
q{dA) belonging to the ground state ensemble of the phase a the statistics of external 
contours in the contour model coincides with that of the true particle model. Note that 
the statistics of non-external contours are different in the contour and particle models. For 
example, the contours of the opposite phase, —a, never appear in the contour model. As 
we explained before the absence of the matching rules makes the analysis of the contour 
model much easier than the analysis of the initial particle model. The price paid for this 
simplification is a rather involved expression for VF(r|g('^ ^^^^). 

s3.3 Statement 3.3 

For all 7 small enough and all chemical potentials X G (A(/3) —7", A(/3) +7"), a> 1/2 
there are /+,a,7 ? ^^sp. f-^\^^, such that for any sequence of cubes A — > IR^ and any sequence 
of boundary conditions q^^''^ belonging to the + (resp. —) ground state ensemble 

^i^, ^log^7^/3,A(Ak-^^^^) = fix,, (3.39) 

Moreover there exists A(7, /3) such that 

3.i0e2 lim A(7,/3) = A(/3) (3.40) 

7— »o 



3.10e3 



and 

The statement is proved in Section 6 using cluster expansion techniques, together with 
a rather explicit representation of the airxiliary pressures f± x-y 

In analogy to the first term Ti in the decomposition (3.25) of the statistical weight 
we introduce the function 

3.7.2 w{r\q^^^^^^^) = (3.42) 

jQ(Supp(r)\5(r)) 'n{<l)=V 



5?^fnK fl,, . (r)) . (^'^(r)) 



g2^TW_^(r)l«_o.(r) ) 



jQ(Supp(r)\5(r)) T?(g)=a(r) i-l 

In Section 4 we will prove the Peierls bound (which is the analogue of (3.30)): 
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s3.4 Statement 3.4 

There is c> such that for all 7 and A e (A(/?) - 7^^, X{P) + 7^^), a > 1/2 

3.7.3 w(r|g('5^(^))) <exp(-c£^£3 V|Supp(r)|) (3.43) 

Finally in Section 6 we will prove 

s3.5 Statement 3.5 

There is c > such that for all sufEciently small 7 and A e (A(/3) — 7", A(/3) + 7°), 
a > 1/2 

logW^^(r|^('5=(r))) - \ogw{T\q^^^^^^^) 

- Yl (/5|I^t^(r)l[/-.(r),A,,-/a(r),Aj)| (3.44) 

m: a^(r)7^a(r) 

<c|Supp(r)|7V4 

Note that the factor 7^/'* is not optimal but it is enough for our purposes. 

By combining the above statements and choosing X — X{P, 7) we then obtain that the 
auxiliary and the true partition functions are equal and consequently that the probability 
of a contour in the original system is bounded as in (3.43). This proves the existence of two 
different phases and together with the exponential decay established in Sections 5 and 6 
concludes the proof of Theorem 2.1. 



3.10e6 



4. Peierls Estimate 



In this section we prove Statement 3.4, i.e. we compare the partition functions in the 
numerator and denominator of (3.42). 

The partition function in the denominator of (3.42) is taken over ground state con- 
figurations of the phase a = (j{T) placed in the region Supp(r) \ 5(r) with the boundary 
condition q^^^^^) specified on 5{T). This boundary condition, different in 5 (F) 

and in 5^{T): it is a configuration of the ground state ensemble of the phase a in the 
region 5^(F) and it is the standard configuration of the phase a in 5^(F). 

The partition function from the numerator of (3.42) is taken over contour configu- 
rations, i.e. ones compatible with if' , in the same region Supp(F) \ S{T) with the same 
boundary condition (F) but with the different boundary condition given on 

5^(F). Contrary to the denominator this boundary condition is the standard configuration 
of the opposite phase, —a. 



Scheme of Estimate 

The estimate of the ratio in (3.42) is performed in several steps. At each step we 
achieve a certain simplification at the price of a non essential error which is a factor not 



18 



exceeding exp ^7"£2^3 '^C^|Supp(r)|^ with a > 0. For 7 small enough the product of the 

finite number of such factors is, of course, dominated by exp (-c^f^g ^C^|Supp(r)|). 

At the first step of the estimate we replace every configuration q = [qi) contributing 
to the numerator or denominator of (3.42) by [gj^^^^^ We show that the error caused by this 
transformation is not important. The calculation is straightforward for the denominator 
of (3.42) since in the contributing configurations the density of particles is bounded. One 
needs to be more careful with the numerator of (3.42). For some of the contributing 
configurations it may happen that in some boxes C(^i) e Supp(r) \ 5{V) the number 
of particles is too large. In this case we replace the configuration inside this box by a 
configuration with a bounded density and we use superstability to check that the error is 
small. As soon as all consideration are reduced to the case of bounded density it becomes 
clear that one needs to check (3.43) for A = A(/?) only. Indeed, varying A in the interval 

A e (A(/3) - 7", A(/3) +7") with a > 1/2 produces an error in Hy^xiq\q^^^^^^\ q±^lr)^) 

which in absolute value does not exceed 7-'^/^|Supp(r)|. Hence the total error in (3.42) is 
negligible with respect to the right hand side of (3.43) 

After shifting particles into the centers of the corresponding boxes integrals over the 

particle configurations q become sums over density configurations g^^'^\q) = {Qx^\q))- 
Accordingly the energy of the configuration q is replaced, again up to non essential error, 
by the discrete version of the mean field free energy functional (3.1) where the function 
Q{r) of a continuous argument r e is replaced by its lattice version Qx^\ x G Z;^^. 

The next observation is that the entropy factor coming from the summation over 
different density configurations g^^'^^ can be neglected and one can consider only a contri- 
bution of two density configurations providing the minimum of the free energy functional 
in the numerator and denominator of (3.42) respectively. This reduces everything to a 
variational problem for the discrete version of the functional (3.1). 

Utilizing (3.16) one can see that for the minimizing density configuration of this vari- 
ational problem the influence of the boundary conditions decays exponentially and at the 
distance of order -y"^"":^ from A*^ it practically disappears. Thus, modulo non essential 
errors, we can replace the initial variational problem by a similar one in a smaller region 
with standard boundary conditions only. 

At the flnal step we consider for the denominator of (3.42) a suitably reduced region 
Supp(r) \ 5(r) with the standard boundary condition Qi3^a- The minimizing density is 
clearly the same gp^a continued inside the region. Investigating then the minimizing density 
configuration for the numerator of (3.42) one takes into account the fact that there are 
sufficiently many boxes C*^^^) inside Supp(r) \ 5{T) where the indicator function I ^ ^ ^ 

forces the minimizing configuration to be different from £»/3,±. This leads to an excess 
of energy, with respect to the energy of Qp ±1 which is enough to produce the factor 
exp(-c£^£3-'^C'|Supp(r)|). 

We now turn to the details. 

Stability Estimates 

First we check that the accumulation of an infinite number of particles in a bounded 
region is impossible. This is a consequence of 
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Lemma 4.1 

For any particle qq and any conGguration q = {qi,q2, . . .) the interaction energy 

Uj{Qo\q) = - ^ J^^Kqo, Qii) + ^ J^'^\Qo,Qii,Qi2iQi3) (4-1) 

is bounded from below 

U^{qo\q) > i^min (4.2) 

by an absolute constant 

i?„,i„= min ^(X^ -S^'^X^ + 2^^'^X) -X (4.3) 

Af>0, 1>7>0 3! ^ ^ 



Proof. Given q let 

^f(r, q) = ^dy^j^ /4_4N 
be a mean number of particles of q situated inside B^{r), r e R^. Then 



- y 4^Kqo,qi,) = -7"" y fdr H ,1 



qii^q qi&q' 

= -^'^ [ dr -a. ^ Mir, q) 

Similarly 



d 



Hence 



(4.5) 



y 4^^(90, 9*3) =7'' / dr I _ 

' J \r-qo\<1 Rd 

qii,qi2,ii3^i (4.6) 

X ^ {-^'{r, q) - S^^M\r, q) + 2^^''N{r, q)) 



Uyiqo\q) I dr H^^ = H^^^, (4.7) 



which proves lemma. ■ 

Lemma 4.1 provides the lower bound 

H^Mq^'^^Iq^'"^^) > ^min|g(^^| = H^,^\A\g, (4.8) 
where g = Igf^^^l/IA] is the corresponding density. To obtain an upper bound take some 

positive £>niax > Qf3,+ - 
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Lemma 4.2 

Consider configurations q^^^ and q^^ ^ such that q^^^ U q^^ -* has at most Qma.x{'^l~^)^ 
particles in any intersecting A cube with the side Then 

H^AQ^'^^IQ^'''^^) < ^^max(^max)|A|, (4.9) 

where 

iymax(^max) = \X{P)\Qm.. + 2^'' Q^,^ (4.10) 



and 



Proof. It is clear that |g(^) | < |A|^max- The strength of four-body interaction between any 
four particles is less than 7^^^. The number of interacting quadruples of particles such that 
one of the elements of the quadruple is a given particle is less than (^'max2'^7~'^)'^ Hence the 
total four-body interaction contributing to Hj^x{q^^^ |q(^')) is less than \A\g^^^2^'^. The 
two-body interaction is negative and does not contribute to the estimate. ■ 

Bad Boxes 

In this subsection we treat boxes containing too many particles. 
Lemma 4.3 

Consider a box C^^\ i < \'y~^ , a configuration q E Q and an integer 

iV > |C(^)|e'^^ = (4.11) 

Then 

dq I ^-(3H,Mi'oW'm < (4.12) 

Proof. We check only (4.12) as it easily implies (4.13). According to Lemma 4.1 the 
interaction [/^(q*-*^' '-'Iq) between particles of q^^^ and q satisfies the estimate 

U^{q^'''''^\q)>H^,^N (4.14) 

It is clear that 

as the maximal volume of the intersection of two balls of radius ^~^Ra is 7"^^ and the 
minimal volume of the intersection of four such balls with the centers in C^^^ is larger than 
^7"^^. Thus the logarithm of the left hand side of (4.12) does not exceed 

- NlogN + N + Nlog\C^^^ - (3H^ir.N 

+ + PY^^ - ,_L,3. iV(A--l)(^^-2)(A^-3) (4.16) 
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It is not hard to check that for sufficiently small 7 and sufficiently large absolute constant 
C(4.ii) the last expression is smaller than —N as soon as A'^ > IC*-^-* |e^('' i^)^. Note that the 
most dangerous term in (4.16) is P'j'^ ^^^^^^ . It is dominated by P^'f^'^ n{n-i){n-2){n-3) 
for > 07-'^ and by iVlogTV - iVlog|C(^)| for |CW je'^f^.iD/S < N < 07-*^. This implies 
(4.12) and hence (4.13). ■ 

Set ^max = 26*^(411)^. An easy consequence of Lemma 4.3 is 

s6.3 Lemma 4.4 

For any contour F and any £ < £2 the partition function in the numerator of (3.42) 
does not exceed 



jQ(Supp(r)\5(r)) vi<l)=V^, Q^'^\q)<Qr, 



m: am{T)i^cr{T) 

6.13 (4.17) 

Proof. To obtain (4.17) we perform a partial integration over the configurations in the 
bad boxes. Suppose that q is fixed outside the cube C'^^^ G Supp(r) where q has more 
than £»max|C'^^^| particles. Now applying Lemma 4.3 we integrate over all such q-s leaving 
exactly ^£»max|C'*'^-'| particles in C^^\ That means that we integrate over the particles in 
excess to |^max|C'^^-'| of them assuming that the total number of particles in C^^^ is larger 

than ^'maxIC''-^-' I- This produces an extra factor 1 + 2e~l'-"^ '''^ < e. We continue this 
procedure box by box and the total number of boxes in Supp(r) is £~'^|Supp(r)|. ■ 

Reduction to Variational Problem 

It was shown in the previous subsection that we may restrict our considerations to the 
configurations with the bounded density. Now we estimate an error which is produced by 
shifting particles in such a configuration into the centers of the corresponding boxes C^^^ . 

s6.4 Lemma 4.5 

Taife £ < 7-^ and configurations q^^^ and q^^"^ with g^^\q^^^ U q^^^^)) < ^^ax- Then 
|if,,,(g(^)|g(^^)) - H,,,{g('\q(^W''^^)\ < ^7|A|(2'^^Lx + 2^'^Lx) (4-18) 



Proof. It is clear that \q^^^\ < lAj^imax- Given two interacting particles the absolute 
value of the error in their interaction due to shifting these particles at the centers of 
the corresponding boxes is less than •y'^i'^. Given four interacting particles the absolute 
value of the error in their interaction due to shifting these particles at the centers of 
the corresponding boxes is less than ■-y^'^£j. As in Lemma 4.2 the number of interacting 
quadruples of particles such that one of the elements of the quadruple is a given particle 
is less than (£»max2''7~'^)^. Similarly the number of pairs of interacting particles such that 
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one of the elements of the pair is a given particle is less than (£»max2'^7 Hence the total 
error is less than ^7|A|(2'^^^^, + 2^'^^^^^. ■ 

This lemma allows us to replace the integrals over dq in the numerator and denomina- 
tor of (3.42) by the sums over density configurations ()^^^\ Namely, consider the partition 
function on the right hand side of (4.17). The integral over q e Q(Supp(r)\5(r)) ^^^j^ 
v{q) — cind Q^^^\q) < Qma.x can be calculated in two steps. 

First one can fix a density configuration g^^^^ — (qx^^), x G [Supp(r) \ d{r)]^^'^^ and 
integrate over configurations qf(Supp(r)\'5(r)) ^^^j^ Q^^^\q) = 

Afterwards one can sum over all q^^'^^ < Qmax compatible with r]^ . The obvious 
upper estimate for this sum is the total number of density configurations Qx < ^max, x G 
[Supp(r) \5(r)](^i) times the maximal contribution given by a single density configuration. 

The number of density configurations g^^'^^ < Pmax compatible with r]^ is less than 

|Supp(r)\<ur)l 

6.19 (^max|C(^^)|) l^^^^'l (4.19) 

Given g^^^'' the integral over configurations q with g^^^^q) = g^^^'' does not exceed 

m: <T^(r)5^<7(r) 

+ log|C(^i)| Yl Qi''^\C^''^\ 

6.17 xe[Supp(r)\5(r)](^i) (4.20) 

y: iog((,(^^)|c(^^)i)!) 

+ /?47|Supp(r) \ 5(r)|(2'^^Lx + 2^'^Lx)) 

as follows from Lemma 4.5. 

Suppose that the density configuration g^^^^ gives the minima of 

2 t^7W-a(r)|y-a(r) ) 

m: o-^(r)#<T(r) 

_|C(^0|iog|c(^i)| Y si''^ ^^''^^^ 

a;G[Supp(r)\5(r)](^i) 

+ E log((^i^^)|C(^^)|)!) 

a;G[Supp(r)\5(r)](^i) 

among all g^^'^^ < £»max compatible with 77^. Then 

+ il' l0g(^^ax|C(^^) I) + /?47(2'^Lx + S'^^^Lx)) 

6.20 (4.22) 
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is the upper bound for the log of the numerator of (3.42). Note that the term 7^/^ in 
(4.22) is the price paid for taking A = A(/?). 

Similarly if ^'^^i) gives the minima of .7-"^^^ ;^(^)(£>(^i^|g('^ ^^^\qi^ ★) among all den- 
sity configurations g^^^^ from the ground state ensemble of the phase a then 

-^.,^,A(/3)(^(^^)k-^'^^"», ?r - |supp(r) \ 5{T) I (7^/^ + /?£i7(2^^Lx + 2''qLj) 

6.21 (4.23) 
is the lower bound for the log of the denominator of (3.42). 

The notations with ~ and ★ foresee forthcoming simplifications and variations. In 

particular -F^,^,A(;3)(^('^)|g(^^(^», gS'^^"^^^) is defined as ^,,pmpM''^\Q^'^^^^^ ^ Q^^'^^'^^^*) 
without the term f E^: .^(r)/.(r) ^^tI?- j(r) k-r/r?^) (4-20). 
Now denoting b-y,i = 1 1 [S^(-)]^^^ 1 1 let 

6.21.1 (xi,X2) = \\[B^{xi)f^ n [B^{x2)f^\\ (4.24) 
and 

6.21.2 /W(a;i,...,a;4) =6-^ ||n|_i [B^{xj)f^\ (4.25) 

be discrete versions of Jj^\xi, x^) and Jj^\xi, . . . ,X4). For any region [A](^) (we do not 
exclude the case [A](^) = Zf) and any density configuration g^^\A) = (qx^), x e [A](^^ 
define a functional 

6.22 E 2^2)^1? (4.26) 

' a;2,a;i€[A]«^) 



6.23.2 



+ ^ E 4^^(a;i,...,a;4)^£^---^i 

a;i,a;2,a;3,K4e[A](^) 

which is a discrete analogue of the mean field free energy functional (3.1). We also define 
a conditional functional 

:F,,p,xiQ^'\A)\Q^'\A^)) = T,,p,x{q^'\A) + ^W(AC)) _ ^^^^ ,(^W(AC)) (4.27) 

with the boundary condition ^(^^ = {Qx^^)-i x G [A^]^^^. Here the similarity with (2.11) is 
obvious and the meaning of q^^^ + q^^^ becomes straightforward if we set q^^^ = 0, a; G [A'^]*^^) 
and = 0, a; e [A]^. Setting 

we introduce 

:F,,0,x{Q^'\A)\g^'\A'^),^) = .^,,^,,(^(^)(A)|^W(A^)) - ^C/,,/3(^^'l(A)|^^'l(A^)), (4.29) 



where g^j^^^ = Qp^^ and u = + or cr = — . With some abuse of notation we use for 
J''y,l3,\{Q''^\^)\Q^^\^^^^^)) an alternative notations T^^i3,x{Q^^\A)\q^^^'^) or 
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■^7,/3,a(^^^HA)|^(^)(5(^^))). The role which is played by functional (4.27) is clarified in 
s6.5 Lemma 4.6 

Consider a V^^^ measurable region A with the boundary condition q^^^^ which is a 
ground state conGguration on every connected component of dA. Then for £ < 7"^ and 
any g^^^ = {qx^), x G [A]^^) such that gi^ < ^max one has 

6.24 (4.30) 

|^7,/3,A(^('nA)|g^"^\*) -^,,/3,A(^(^nA)|^('n9^^^^),*)l < ^75(2'^^Lx + 2^'^^Lx + ^Lx)|A| 

6.24.1 (4.31) 

Proof. Estimate (4.31) is an obvious consequence of (4.30) and we concentrate on (4.30). 
The difference between ^7,/3,a and .^7,/3,a has two sources. The first one is due to re- 
placement of balls -B'y(-) in the definition of jj'-* by their lattice versions [B^{-)]^^^ in the 
definition of Clearly the difference between |-B7(-)| — and | [!?-),(■)] *^^)| is less than 
7~'^£7. Hence the error produced by the discretization of -B^(-) can be estimated exactly 
as in Lemma 4.5. 

The second source is due to not properly counted contribution of pairs xi,X2 with 
xi = X2 and quadruples xi,X2,X3,X4 with not all Xi^ being different. To estimate from 
above the absolute value of this error let us consider the following five contributions to the 
energy of g^^\ 

(i) Self-interaction of ci^^ due to the pair interaction of particles in ci^^ . 

(ii) Self-interaction of Cx^ due to the four-body interaction of particles in cjif^ . 

(iii) Interaction between Cx^ and Cx^ due to the four-body interaction of two particles in 
Cx^ with two particles in Cx^ ■ 

(iv) Interaction between ci^^ and ci^2 due to the four-body interaction of three particles 
m Ci^^ with one particle in ci^^ ■ 

(v) Interaction between Cxi , CxJ and Cx^ due to the four-body interaction of two par- 
ticles in cjif^ with one particle in ci^^ and one particle in ci^^ . 

All five contributions can be estimated by similar arguments. For that reason we present 
these arguments only for cases (i) and (v). 

The strength of self-interaction of Ci due to the pair interaction of particles in Ci 
is less than 7'^(Pmax^'^)^- The number of boxes Cx^ in the region A is £~'^|A|. Hence the 
total contribution is less than {i'y)'^g'^g^-^\A\. 

(£) (£) (£) 

The strength of interaction between Ci/, Cx2 and Cig due to the four-body interac- 
tion of two particles in Cx^ with one particle in Cx^ and one particle in ci^^ is less than 
7^^(^max^^)^- The number of boxes C^t^ and Cxs interacting with given box Cx^ is less 
than {'y£)~^'^. The number of boxes ci^^ is £~'^\A\. Hence the total contribution is less 
than (£7)'^max|A|. ■ 
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Looking for the minima of ^^■y,i3,\{Q^^^\Q^^\q^^^'')) it is simpler to understand gi' as 
continuous variables. Therefore if the minima of J^-y.f3,\{Q^^^\Q^^\q^^^^)) is achieved on 
density configuration g^^^ then it may happen that at least for some x G [A]^^^ the number 

i'^gx^ is not an integer. The solution to this problem is given by 
Lemma 4.7 

For the density conhguration gi^ = £~'^[£'^gx^] one has 

|^7,/3,a(^(')(A)|^(^)(q(^^))) - .^^,^,a(^(')(A)|^(^)(q(^^)))| < ce-^grr^.^lM (4.32) 



(Here [ • ] denotes the integer part of a number.) 

Proof. We proceed as in the proof of Lemma 4.5. The total number of points x G [A]^^) 
is £~^\A\. Given point x the absolute value of the difference in the corresponding self- 
interaction is less than I3~^ log ^max + |A|^max- The difi^erence in two-point interaction be- 
tween points xi,X2 G [A]*^^-* in absolute value does not exceed 27'^^max^'^- The number of 
points X2 interacting with given xi is less than (7^)"'^. The difi'erence in four-point interac- 
tion between points xi, X2, X4, X4 G [A](^) in absolute value does not exceed 47'^'^(^max^'^)"^- 
The number of points X2,X3,X4 interacting with given xi is less than {■j£)~^'^. Combining 
this estimates one obtains the lemma. ■ 



Variational Problem (Dependence on Boundary Condition) 

The results of the previous subsection reduce the Peierls estimate to the following 
variational problem: 

(i) Find the minimum of T^^f3^x(p){g^^^\Supp{T)\5{T))\q^^^^^^\ qi^^^^^\ic) over the den- 
sity configurations g^^'^^ = {g^x^^)-, x G [Supp(r) \ 5{V)]^^'^^ such that g^^'^^ belongs to 
the ground state ensemble of the phase a. 

(ii) Find the minimum of J^^,,^j,;^(^)(^(^i)(Supp(r) \5(r))|g(^~(^)\ gi^^ ^^^\^) over the con- 
figurations = {,Qx^^)i X G [Supp(r) \5(r)]*^^i) such that g^^^'' is compatible with 

(iii) Estimate from below the difference between the minimal value of 

^T,/3,A(/3)(^^^^HSupp(r) \ 5(r))|g(^=(n), gi^^(^»,^) and the minimal value of 

^7,/3,A(/3)(^(^^Hsupp(r)\5(r))|g(^^(n),gr('^»,*). 

The existence of the minima above is obvious as < ^li; < £»max- Following the approach 
of Section 3 define 

n,,£{x,g('^) = b-} Yl (4-33) 

xie[B^ix)]W 



and 

(e) 

I,4x,g('^) = b-} Y: ^(logp£^-l), (4.34) 



xt&[B^{x)]W 
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where we deliberately use a general scale {£) instead of (£i) as the constructions below are 
of general origin. In complete similarity with (3.4) 

6.28 (4.35) 
implying that g^^\Z'^) = Qp,x,- ^ind Q^^^'Lf) = Qf3,x,+ are the global minimizers for 
^■y,/3,x{Q^^^)- The local minimizers of J^'y,f3,x{-)i i-e. the minima of J^j^j3^\{-\g^^\A'')) or 
^'y,p,xi'\Q^^\-^^):'^): cire studicd in the lemma below. Note that in this lemma we consider 
not only A — A(/3) but all A G {X{/3) — S{/3), X{/3) + S{/3)). Though the lemma discusses 
^'y,p,xi-\Q^^HA^)) the same argument covers the case of J-'j^i3^x{-\g^^\A'^),'k). 

s6.7 Lemma 4.8 

Consider a V^^^ measurable region A and take a boundary condition g^^\A^) with 
max I g^^^ — gs x a\ < where a is one of the phases, + or —. Then the unique minimum 

xe[dA]W ' 

of T^^I3^x{q^^\-^)\q^^\-^'^)) ^^ong density conhgurations g^^\A) with max \gl^^ — g)3,x,a\ ^ 

a;e[A](^) 

( is achieved on the density conGguration g^^^ = {gi^), x e [A]^^^ such that 

l^i'^ - QP,xA < ^^^^5^^:,^ , (4.36) 



1 - a{(3) 



Proof. Calculating -^^i^J^-y^p^x{g^^\A)\g^^\A'^)), x e [A](^) one obtains the necessary 



where [•] denotes the integer part. 
Proof. Calculat 
condition for ^^^^ 



0=^log^(^)-A-j;/f (x,xO^£) 

6.30 (4.37) 

+ ^ E 4'n-,-i,-2,x3)^i?^S^S, 

It is clear that for gi'' = gp^x^cn x G [A'^]'^^-' one has gi^ = ^/3,a,o-, x G [A](^\ Introduce 
an aiixiliary parameter t G [0, 1] and an interpolated boundary condition 

6.31 gi'\t) = (1 - t)gp,x,a + tgi'^ ,x G [A'^]^^) (4.38) 
Let g^^\t) be the solution of (4.37) with the boundary condition g^^\t). Then 



^(f) {t) = exp i(3\ + (3 (^) 

\ a;ie[B^(a;)](^) 



6.32 ^ (4.39) 

I I\'\x,xux„ xs)gi'l it) (t) (0 

Xl,X2,X3e[Bj{x)](^^ 
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for all X e [A] . Taking the derivative with respect to t one obtains 



d 
di 



\xie[B^(x)lW 



[B^{x)]i 

i 4^) (X, xi, a.., X3) (t) 

a;i,a:2,a;3G[S-y(a;)](^) 



(4.40) 



Denote 

I^^\x^,X2\Q^'\t)) = I^'\xi,X2,x^,x^)e^}{t)e^l{t) (4.41) 

It is clear that 

(xi,X2)(^;3,A,a-C)' < (xi , X2 1 ^(^H^)) < (xi , X2) (^;3,A,a + C)' (4-42) 

as soon as max \ Q^^\t) — Qi3,x,a\ < C- 
Introduce a symmetric matrix 

A{x^,X2\Q^'\t)) = -I^^\x„X2) + lI!,^\x^,X2\g^'\t)) (4.43) 

and a diagonal matrix 

D{x^,X2\g^'\t)) = {Pgi'l{t))-'TL , (4.44) 

Xl=X2 

where Xi,X2 G [A]*^^^. It is not hard to see that the inverse matrix B = {D — A)~^ exists 
if max \gl^\t) — gf3,x,a\ < C- Indeed, for such gx\t) this matrix is given by a convergent 

series 



\n=0 



B = ( J](L>-iyl)'= ) D-^ (4.45) 

because 

\\D-'A\\=m^ ( Pgi'^{t)J2\M^u^2\g^'\m 



(4.46) 

<a{p) 

<1, 

where a(/3) is defined in (3.15). Moreover, the representation (4.45) and the fact that 
A(xi, X2\g^^\t)) = if dist {xi,X2) > imply that 

\B{xi,X2\g^'\t))\ < ^(^^^A:^±0a(/3)[7dist (4^47) 
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Iterating (4.40) and observing that 4fQx\t) = —Qi3,x,a + Q^x^ for x e [A'^]^^) we rewrite 



(4.40) as 

d 
dt 



E E Bix,x^\g^'\t))Aix,,X2\g^'Ht)M^-Qf3,x,a) (4.48) 



a;ie[A]W a:2e[A=]W 

The right hand side of (4.48) is an absolutely convergent series in terms of Qx\t)-, x G [A]*^^) 
and solving this differential equation one finds Qx\t)- The solution exists at least up to 
t = t\ at which the condition 

d 



max 

x&[A]W 



dt 



Q^x\t)-Qf3,) 



is violated. Suppose that ti < 1, i.e. 



max 

x€[A]W 



X.a 



for t < ti and 



d 
dt 



Q^x}{tl)- Qp,X,cr 



C 



for some Xi G [A]*^^-*. Then the representation (4.48) is valid for t G [0, ti] and 



max 

x&[K]W 



for t <ti as follows from the obvious identity 

ft ^ 

Plugging (4.50)-(4.52) into (4.40) and using (3.16) one concludes that 

d 



max 

x£[A]W 



dt 



Qi^\t) - Qf3,X,a 



for t <ti which contradicts (4.51). Hence 



max 
xe[A]W 



and 



max 

x&[K]W 



(4.49) 

(4.50) 
(4.51) 

(4.52) 

(4.53) 

(4.54) 

(4.55) 
(4.56) 



for all t G [0, 1]. Thus representation (4.48) and estimate (4.47) are valid and being joined 
with (4.53) they give (4.36). 

The density configuration ^(A) was defined as the solution of (4.37) and we just 
checked that such a solution exists and satisfies (4.36). On the other hand, it follows 
from (4.42) and (3.16) that the Hessian matrix ■^-j^^j^J^^^f3^x{g^'^\A)\g^^\A'')), x,y G 

[A]*^^) is positive for any ^'^^^(A) G (^»/3,a,o- — Qi3,x,a + C)^ with the mass bounded from 
independently on ' (A) and ^(^)(A^). Thus the function J^^^p^x{-\Q^^\^'')) with convex 
domain {g/3^x,a ~ C; Qi3,x,a + C)^ is convex and therefore ^(A) is its unique minima. ■ 
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Variational Problem (Comparison of Two Minima) 

In this subsection we return back to the case A = A(/5) and £ = £i. First we need a 
stronger version of Lemma 4.8. 

S6.7.1 Lemma 4.9 

Consider a measurable region A and take a boundary condition g^^'^^A^) with 

max l^lf^-* — Pa X al < C? where a is one of the phases, + or — . Then the unique mini- 

mum of J^^^I3^x{g^^^^A)\g^^^\A'^)) among density conGgurations g^^^\A) with max \gi^'^^ — 
Q(3,x,a\ < C achieved on a density configuration g^^'^^ — (gx^^), x e [A]*^^i^ such that 

- QP,x,a\ < ^ , (4.57) 

where [■] denotes the integer part. 

Proof. The difference from Lemma 4.8 is that we require here closeness to gp^x,^ only 
on the scale £2-, which is much larger than £\. As a consequence the variables g^x^^ may 
have larger fluctuations and the free energy functional ^7,/?, a is no longer convex. The last 
argument in the proof of Lemma 4.8 is then not valid anymore. We will prove that any 
minimizer gi^^ (we do not have yet uniqueness) satisfies the bound 

6e.29c l^i^'^ - ^/3,A,a| < C (4-58) 

After (4.58) the proof becomes essentially the same as for Lemma 4.8 and it will be omitted, 
we will just prove (4.58). To this end, it is enough to show that if ^if^-* > g/3,x,a + C) 
X e [A](^i), then 

-^^,,^,x{Q^''Hm'^''\^')) > (4.59) 
(ygx 

and that the reverse inequality holds if g^^^ < g(3,x,a ~ C- Since the two proofs are similar 
we only consider the former. As in (4.37) we have 



6e.30 



6e.31 



d 



dgi 

where 



^J'j,p,x{Q^'^Hm^''\^n) = ^logg^J^^-X-G^{x,g^'^^) (4.60) 



G7(a:,p(^^)) = J]/^^,\(a:,a:i)p(^^)-i J2 I^%{x,xi,X2,xs)gi\^^ g^^,^^ gi'i^ (4.61) 

Xl ' Xl,X2,X3 

with the I^£^ as in (4.24) and (4.25), having added the subscript £i to specify the scale 
used in the definition. Since we are supposing that gx^^ > ^/3,a,ct + C; there is c > so that 

loggx'^ > log(g^,A,a + C) > logg/?,A,a ^ C _ ^^2 

/3 ~ P ~ P PQ(i,X,a 
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We next prove that there is c' > so that 



3 2 

G,{x,Q^'^^)-[qp^,.-^-^] <C|1-^I + c'(7'^^ + C') (4.62) 



2 



This is obtained by adding and subtracting I^l^ to I^}^ and then expanding. In the terms 

containing I'f^} we can partial-sum the densities Qx^^ over x in the same cube C^^'^\ as 

I'^l^ are constant in such cubes. We thus reconstruct a density q^^'^^ which, by hypothesis, 
differs from gp,\^cr at most by C. In this way we obtain the square bracket term in (4.62) 
and the terms with on the r.h.s. of the same equation. The last term on the r.h.s. of 
(4.62) comes from the the difference between the /^"^^'s. We have thus proved (4.62). 
By recalling that Qf3,\,a satisfies the mean field equation, we get 

-^.F,,^,,(^(^^)(A)|^(^^)(A^)) > — ^ - C|l - ^1 - cC' - (4.63) 

Since the coefficient that multiplies C on the r.h.s. of (4.63) is strictly positive (by the 
choice of /? and A), the r.h.s. of (4.63) is positive for 7 and ( small enough. (4.59) is thus 
proved and since the proof of the reversed inequality for g^^^ < Qi3,\a ~ C is completely 
similar, we can then conclude that \Qx^^ — Qp,\,a\ < C. The proof of the lemma is now 
analogous to that of Lemma 4.8 and it is omitted. ■ 

For any contour F the set 

S = a(^3)Ext(r) U j U d^^^hnimiT) I (4.64) 

\m: a^(r)=a(r) / 

is occupied by the ground state of the phase cr = (t{T). Consider a strip 5^(r) of the width 
situated in the middle of S. According to Lemma 4.8 inside 5^(r) the density con- 
figurations minimizing JF^^^^a(^»*-^^'' 1?*-'^ ^^^\q^ 8jaA T^^p^\{Q^^^^\cf^ ^^^\qi.„^^^\i<) 
differ from Qj^^^ by at most (1 — a{l3))~^ C,a{l3)^'^'^'^ Therefore the density configura- 
tion which minimizes J^'y,i3,\{Q^^^''\(f^ ^^^\q^ *"'"'*'',*) under the additional condition that 
this configuration coincides with Qfs^^ in 5^(r) gives the value of the true minima of 

^'y,l3,x{Q^^^^\Q^^ '■^■'^ '■^■'''i*) up to nonessential error. The same is true for the density 

configuration minimizing J^'y,j3,\{Q''^^^\q^^ Q'^o- 5 *) with the same additional condi- 
tion. Obviously these modified minimizing configurations coincide with each other not 
only in 5^(F), where they both are equal to gp^aj but in the whole part of S stretching 
from 6={T) to S={T), i.e. in 

= 5(W2)Ext(F) y I y d^^^/'^hntmir)\[j5={T) (4.65) 

\m: a^(r)=a(r) / 
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6.52 



Estimating from below the difference between the minimal value of ^7,/3,A(^^^^^k^^ ^^^\ 

q^^cr ^^^\'^) ^iid the minimal value of J^^,i3,\{Q^^'^^\q^^ Qa ^^^\*) the contributions cor- 
responding to Si cancel each other. Hence the initial variational problem in Supp(r) \5(r) 
with a general boundary condition imposed on S{r) is reduced to a similar problem in a 
smaller volume Supp(r) \ Si. 

For the sake of notational simplicity from now on we suppose that the boundary 
condition is a standard one already for the initial problem, i.e. it is equal to gp^^ in 5^(r). 
Thus all we need to finish the Peierls estimate is the lower bound for 



a;e[Supp(r)\5(r)l(«i) 



6.49 (4.66) 

Wc note that for x such that dist (x, 5{T)) < 7"-^ the values of Xy,£(.T, g^^^) and TZ^j{x, g^^^) 
depend on the boundary condition which is equal to gp^a in ^^(r) and in 5^(r). 

s6.8 Lemma 4.10 

For any contour T = (Supp(r), r]^) 

6.50 min AJ^^,p,^{g^'^\T) > ceie^''e\Snpp{T)\ (4.67) 
Proof. Rewrite A.7^y,^,A(^'*■^^^ T) in the form 

a;e[Supp(r)\5(r)](^i) V 

+ E ^^oggiJl^ (4.68) 

xie[B^,(x)]{^i) 

_ThAi^f^logn,M^,g'^'^^) 
and observe that 

6.51.1 Fp,x{n^,iAx, g^^'^)) - F0,x{Q0,a) > (4.69) 
by definition of gf^^^ and 

E ^log.g^) > ^^'^-("'^^'^^^ og7^,,,,(x,,(^^)) (4.70) 



by convexity. 
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We call Ci^^^ e Supp(r) a wrong box if either 77^ = or rj^r]^^ = —1 for at least one 
cube Cx^^ adjacent to ci^^\ According to the definition of the contour there exist at least 
(3-^3)~'^|Supp(r)| wrong boxes cjif^^ e Supp(r) such that they are at the distance longer 
than 57"^ from each other. In particular C^^^ ^ fl C^^'^ ^ = for any two such boxes 
Cxi ^ and Cx2 ^ • 

Consider a wrong box Cx^^ for which rj^ = 0. If inside Ci^'^ there exist at least 
^r^^2 points xi with 

n^M^uQ^''^) ^ - 1, 0p,- + U - 1, ^^,+ + 1) (4.71) 

then ci^'' ^ contributes to AJ^y^j3^x{g^^^\T) by at least 

die min F^',a(^/3,+)) (4.72) 
This contribution comes from the terms 

eiFf3,x{n^,eAxuQ^''^))-ifFp,x{Qi3,a) (4.73) 

in (4.68). 

In the opposite situation when 

n^,iA^^,g^'^^)^(^gp,.-^, _ + |J (^^^^^ _ ^, gp,+ + ^^ (4.74) 

for not more than i^'^li points xi G ci^'^ ^ we extract a contribution similar to (4.72) 
from the terms 

E - iogR,,.(xi, (4.75) 

x2elB^ix^)](^i) 

in (4.68). 

Consider xi e C^^''' such that Cx'^^ e [5^(a;i)]^^i) and observe that 

,(..) = ||[^C..)](.0||-i Y: (4-76) 

does not belong to (^^ _ - C, ^/j - + C) U (^/3,+ - C, ^/3,+ + C) as r]l = 0. Denote 

Tlx, = \\[B,{x^)f^^ \ [Ci'^^n-' Yl (4-77) 

cc2G[B^(a;i)](^i)\[ci'2)](£i) 
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Then by convexity 

X2e[i?^(xi)](^i) (4.78) 
implying the lower bound 

for (4.75). Now we apply the inequality 

(1 — a)a log a — ah log h— ((1 — a)a + ah) log ((1 — a)a + ah) > — — — — 

^ ' ^ '2 max(a, h) a 

(4.80) 

which is true for any a, 6 > and < a < 1. This leads to the lover bound 



pd 1 1 J 1 1 \ / _ pa [ iiL""^- J II I \ /_ (A o-i\ 

for (4.79). Observing that (&7,^i)"^||[Cx^'^]^^'^|| has the order 7^*2^ we conclude that (4.81) 
exceeds 



(^7^ (^^,,(^0) _ ,(^.))^ > cifr^'^-^ (4.82) 

POj,i^ Qp- \ / PQ(3- 

.(27-') 



for 7 sufficiently small. The number of points Xi € Ci for which (4.82) is true is not 
less than h^^i^/2. Remind that for these points 



(4.83) 



and [5^(a;i)](^i) 9 Ci^^^ Thus (4.82) again gives us a lower bound 

c£iQ\PQp,-)-^ (4.84) 

similar to (4.72). 

Finally for the wrong box C^^'' of the second type, i.e. when ry^ry^^ = — 1 for an adja- 
cent box C^^\ we can consider C^'^'' UClc^^ instead of C^^^ and repeat all the arguments 
above. They will work perfectly because 
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is again outside (^^3 _ - C, Q/s- + C) U (^/3,+ - C, 0f3,+ + 0- ■ 

To finish the proof of the Peierls estimate we need to compare all the errors which are at 
most of the order c£i7|Supp(r)| = 07^1 |Supp(r) | (see (4.22) and (4.23)) with the contribu- 
tion coming from Lemma 4.10. The last is of order c^^^g '^|Supp(r)| = c7^"2+da3|Supp(r)| 
and dominates 07^*1 |Supp(r)| since da2 + da^ < ai. 



5. Auxiliary Model 

In this section we study metastable models and we prove that the corresponding 
measures satisfy the Dobrushin uniqueness condition and hence exhibit an exponential 
decay of correlations. We note again that the configuration q^-^^ in a P^^^) measurable 
region A belongs to the ground state ensemble of the phase a iff in every cube Ci^'^ e A 
the density of particles 

(5.1) 

belongs to the interval 

4-2 {Q(3,a - C, Qf5,a + C) (5-2) 

The partition function of the auxiliary model is given by (3.36), i.e. it is the integral over 
ground state configurations of the corresponding contour partition function. 

In the previous section we have shown that some partition functions initially defined 
in terms of particle configurations can be approximated by partition functions written 
in terms of density configurations related to some scale £. Now we go further and show 
that the auxiliary model for each of two phases can be equivalently rewritten in terms of 
density configurations. Such an equivalent model is defined on the lattice Z^^ with the 

density variables g^'^^ taking discrete values n£^^, n = 1, 2, . . . from the bounded interval 
(5.2). The corresponding Hamiltonian is of infinite range but with sufficiently fast decaying 
interactions. Quite naturally this Hamiltonian is close to (4.27) (with £ = £2) and it can be 
understood as a small perturbation of a positive definite quadratic form. The treatment 
of equivalent model is based on a specific approach [COPP] to the Dobrushin uniqueness 
theorem developed initially for unbounded lattice spin systems. 

Prom now on we fix £2 as the scale at which we define density configurations. Thus all 
regions are assumed to be X)^^^) measurable and we often drop the superscript {£2) from 
notations. 

In the next subsection we construct the effective Hamiltonian the exact form of which 
is stated in Lemma 5.2 at the very end of subsection. Then in the last two subsections we 
prove Lemma 5.3 saying that the effective Hamiltonian satisfies the Dobrushin uniqueness 
condition. 



4.1 



gi'^Kq^'^^) = £2' g(^)nc^^) 
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Reduction to Density Model. 

The technical part of the reduction is based on the cluster or polymer expansion 
technique. For the convenience of the reader Section 7 quotes a version of the general 
cluster expansion theorem which is suitable for our purposes. 

Consider a region A with the boundary condition q^^ ^ belonging to the ground state 
ensemble of the phase a. For the partition function Z^^ ^^{A\q^^ ^) we decompose the 
integral in (3.36) into a sum of integrals. In this decomposition the external sum is taken 
over all density configurations Qx, x G [A] (^2) satisfying (5.2). Given such a density con- 
figuration p(A) an internal integral is taken over all particle configurations q^'^^ such that 
q'x^\(i'^^^) = q'x^\^) for all x G [A]*^^^)^ xhe reduction we perform is nothing but the 
calculation of 



log 



/ dq^^^S.,,, e-^^-^(«''^l^"^'^') y Ww^{rf^\q^'^^^^^^) (5.3) 



as a function of ^x(A). Denote by q^^i the particles of g*^^^ situated inside Cx^"^ and set 
Tlx = 4Qx = q^^^ n dj'^ . Then /g(^) dq is 



n Sf / ■ • • / n ft 

where the integral can be understood as an expectation with respect to the system of 
independent particles qx,i uniformly distributed in the corresponding boxes Cx'^\ 

The Hamiltonian H^,x{q^^^q^^^'^) can be decomposed into the sum of a ^-dependent 
Hamiltonian 

H,,x{Q^'^\<1^''W'^^) = H.^xiiq^^^n^"^'^) = if^,A(^(A)|g(^^)) (5.5) 
and error terms 

Aiy,,,(g(^)|g-(^^)) = H,,,{q(^^\q(^'^) - H,,,{g('^\q(^W''^^) (5-6) 

Clearly exp (— /9iJ^^x(^|9^^'^^)) can be taken outside the integral in (5.4) leaving us with 
the calculation of the log of the partition function 



{rj'^eA i 
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where the first product is included for the convenience of treating this partition function 
as an expectation over a system of independent particles. 

Observe that the error part of the Hamiltonian is given by 

4.8 AH,,x{q^^^\q(^^^) = - J2 Yl AjW(Q,„fe,Q,3,Q,J (5.8) 

where Aj\'^\qij^,qi^) and Aj!f\qi^,qi^,qi^,qi^) are much smaller than 1. The estimates 

and 

4-10 < Aj!f\qi^,qi^,qi^,qiJ < c^''^J^'^\qi^,qi^,qi^,qiJ (5.10) 

are true unless some of the interacting particles are at the distance larger than ^-i— ^-i+"2 
from each other. In the last case AJ^^\qi^,qi^) and AJj^\qi^,qi^,qi^,qi^) are extremely 
small 

< 7-'Ajf (g,,,g,J, ^-'''Aj!f\q,„q,,,q,„q,J < c^^d/2-a2/2 

Denoting 



4.14.1 



4.12 w\^\qi„qi,) = e^^-^^ '^''^.5^2) _ 1 (5.12) 



and 

4.13 w^^^\qi„qi„qi„qu) = g-^^-^^'^^^n .5«2.5»3'<?m) _ 1 (5.13) 

we have 

AAA (5.14) 

First we obtain the polymer expansion for the log of the partition function 

^a;e[A](^2) / \a;e[A](^2) i=l ' / 

containing no contours. Opening all brackets in (5.14) we rearrange the expression under 
the integral in (5.15) in the following way. 

Let 2-link, L^^) = (q'i,q'2), be a couple of particles q'i,q'2 such that W{L^'^^) = 
w'^\q\,q2) 7^ 0. Similarly a L^^^ = {qi, q2, Qs, Qa), is a quadruple of particles 

<?i, 92, q'4 such that VF(L(^)) = w\^\qij q2j qs, qA) 7^ 0. The quantities W{L^'^^) and 
W{L^^^) are called the statistical weights of the 2-link and 4-link respectively. Two links 
are connected if they have a common particle. We stress that if two links have no common 
particles but have a common space point occupied by one particle from the first link and 
another particle from the second link then these links are not connected. A pre-diagram, 
9, is a connected set of links. Denote by q{9) = {qi{6)) the particles of q^^^ which influence 
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9, i.e. the endpoints of links of 9. The statistical weight, w{9), of the pre-diagram is the 
product of the statistical weights of the contributing links. Two pre-diagrams are compat- 
ible if they arc not connected. Finally, a compatible collection of pre-diagrams consists of 
mutually compatible pre-diagrams. 

The definitions above justify the representation for (5.15) of the form 



^a;e[A]«^2) i=l / {6ij}0A= j 

where the sum goes over compatible collections, {9j}, of pre-diagrams and {9j} ^ means 
that every 9j has at least one particle inside A. 

We call two pre-diagrams equivalent if they can be transformed one into another by 
shifting some particles such that every shifted particle q^^^ ^ remains in its initial box Ci^^\ 
The corresponding equivalence classes, O, are called diagrams. To have a geometrical 
interpretation of the diagram we identify it with the pre-diagram 9 & Q having all 
particles at the centers of the corresponding cubes Cx'^K We say that two diagrams ©i 
and ©2 are compatible if any 9i e ©i and 92 G ©2 are compatible. Setting 

^(©) =^-'^l''(®)l / w{9) (5.17) 
Jeee 

we rewrite the partition function (5.15) in so called cluster form 

where the sum is extended to all compatible collection of diagrams. The transition from 
(5.16) to (5.18) relies on the fact that for compatible 9j-^ and 9j^ the corresponding sets of 
particles q(9j-^) and qifij^) do not intersect each other. 

Lemma 5.1 

Let 

a(©) = |Q(©)| (5.19) 

Then 

E \W{Q')\e<'^'^ < a(©), (5.20) 

where Q' Q denotes a diagram ©' not compatible with a given diagram ©. 

Proof. We say that a diagram is not compatible with a particle q and we denote it 
7^ g if g G g(0). Prom our definition of compatibility of diagrams it is clear that (5.20) 
follows from 

E |W^(6>)|e"^®^ < 1 (5-21) 
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Remind that for each box C^'^\ x e [A (^2)] number of particles inside C^'^^ is fixed 
as the density configuration q is fixed. Since q belongs to interval (5.2) there are at most 
c{qi3,\,cf + C)l~'^ 2-links and at most c(p^ > o- + C)'^7~'^'^ 4-links passing through any given 
particle. The statistical weights of links satisfy 

|W^(L(2))| < c-f'^^r^d (5_22) 

and 

|Ti"(L(4)) I < 07^^73^ (5.23) 

as follows from (5.9)-(5.11) and (2.3)-(2.4). 

We provide a diagram with an abstract tree structure according to the following 
algorithm. The root of the tree is the particle q. Links which start from q are called links 
of the first level. Links which start at endpoints of the links of the first level and are 
different from them are the links of the second level. Generally, links which starts at the 
endpoints of n-th level links and are different from all links of levels 1, 2, . . . , n are called 
links of level n + 1 . 

Denote by n(©) the maximal level of links in O. It is clear that for 7 small enough 
|W^(e)|e'^(®) < JJ (l+e2|T^(L(2))|) J]^ (l+e^|W^(LW)|^ _1 < 1 (5.24) 

e/q: n(G)=l L(2)9q L^Bq 

By induction suppose that 

|W^(e)|e«(®) < 1 (5.25) 

Q'/'q: n{e)<N 

and consider © with n(0) < + 1. Take a link of the first level in such 0. From every 
non root endpoint of this link "grows" a subdiagram Oi with n(0i) < N. Hence 

J2 |W^(e)|e"(®) < JJ (l + \W{L^^^)\{e + lf^ 

e/q: n(e)<JV+l L(2)9g 

X II (^i + |H-(LW)|(e + l)4j -1 (5.26) 

< 1 

Here e correspond to the case when nothing is "growing" from a given endpoint of the 
first level link while 1 is the inductive estimate for the case when nonempty subdiagram 
©1 with n(6i) < N is "growing" from this endpoint. ■ 

From Lemma 5.1 applying Theorem 7.1 one obtains the polymer expansion 

E ^(^) (5.27) 
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for the log of (5.18). The precise definition of the polymer tt = [0^^] and its statistical 
weight VF(7r) can be found in Section 7. Geometrically a polymer is again a diagram-like 
object probably with some links entering it more than once. We underline that constructing 
pre-diagrams, diagrams and polymers one considers all particles entering these object as 
distinct, say having unique indices or labels. For that reason we call polymers from (5.27) 
labeled polymers. To clarify the dependence of (5.27) on we perform another factorization 
and define unlabeled polymers. 

Suppose that from the total = Qx^t particles situated inside a box ci^^'* exactly 
k{Ti) particles contribute to the labeled polymer tt. Replacing these A;(7r) particles with 

another A;(7r) particles from the same box ci^^'* one obtains different labeled polymer with 
the same statistical weight. Two labeled polymers which can be transformed one into 
another after several replacements, possibly taking place in different boxes, are called 
equivalent. The corresponding equivalence classes are called unlabeled polymers and are 
denoted r. In other words, unlabeled polymer is obtained from a labeled one by dropping 
out labels of particles. 

Denote by X{t) C [A]*^^^) gg^ Qf ^j^g centers of all boxes Cx'^'' containing particles 
from r. For x e X{r) let kx^r) be the number of particles from C^"^^ contributing to r. 
Then the total number of different labeled polymers tt e r is a polynomial function of 

QJi X e x{t) 

4.27 < P(t) < {Qx^if''^^^ (5.28) 

Setting W{t) = W{7r), where tt is an arbitrary labeled polymer from r, we obtain the 
expression for the log of the partition function (5.15) 

4.28 J2 W{t)P{t) (5.29) 

tGA 

written in terms of Qx- Despite its involved structure we need only few simple estimates 
on this sum. 

As follows from Corollary 7.2 in Section 7 the sum of statistical weights of all labeled 
(and hence unlabeled) polymers passing through given particle and containing not less than 
k links does not exceed ^-'="2/2 jjence the sum of the statistical weights of all labeled 
(or unlabeled) polymers containing two given particles qi and q2 on the distance r > 
from each other 

4.29 ^ \W{7t)\^ Y1 |W^(t)|P(t) <7-1^"]"2/2 ^5_3q) 

71-991, 52 TBqi,q2 

Similarly for any given particle q and sufficiently large absolute constant c 

E I^WH E I^(t)|P(t)<7'', (5.31) 

nBq: L{t)>c rSq: L{t)>c 

where L{t) denotes the number of links contributing to tt or r. 
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Polymer sum (5.29) can be viewed as a Hamiltonian which we separate in two parts. 
The first one is 

4.31 Aiy(;i(^|g(^^)) = Yl W{t)P{t) (5.32) 

r0A<=: L(t)<c 

and 

4.32 AH^%\q(^'^)= Yl W{t)P{t) (5.33) 

T0A=: L(t)>c 

It is not hard to see that AH^j^{g\q^^ )) is simply a finite radius Hamiltonian of the 
polynomial type 

4.33 AJfW (^|^(^^)) =4Y1 ^(^) n (5-34) 

Here the sum is taken over connected sets of links (with no restriction for a given link to 
enter this set more than once) containing less than 0(5.31) links and the product is over all 
endpoints of the links. The notation Qq instead of Qx is not ambiguous as all endpoints of 

the links are assumed to be at the centers of the corresponding boxes C'^^\ The statistical 
weights W{D) are obtained by resummation from (5.32). 

From (5.15) we pass to partition function (5.7) containing contours. For the log of the 
ratio between (5.7) and (5.15) we also obtain a polymer expansion exploiting the theory 
of contour models with interaction [DS] , [BKL] . This expansion has the form 

4.34 Y ^(0 (5-35) 

where ^ are another polymers constructed from labeled polymers tt and contours F in the 
same way as polymers tt are constructed from diagrams G. The statistical weights W^(^) 
are local functions of Qx and we interpret the whole sum as the Hamiltonian 

AH^^l{e\q^^')) = J2 WiO (5.36) 

The only property of this Hamiltonian used later is the estimate 

4.36 E 1^(01 <7''^ (5.37) 

Since all details can be found in [DS] and [BKL] we give only a sketch of the proofs pointing 
out few technically important moments. 

Denote by ^'(•|g^^'^-') the Gibbs distribution of Ylx€[A](^2) '^2Qx particles given by the 
Hamiltonian AHj^x{q^^'^\q^^ -') with the corresponding partition function (5.15). For every 
contour F consider region R{T) = 9*^^^/'^-'Supp(F) with empty boundary condition 0(^('") ) 
and define a modified statistical weight 

4.37 iy^(r/^|?(^=(r))) = i.(5(^=(r))|0(«(r)^)) W^iv^\q^'=^^^^) (5.38) 
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Then the ratio of partition functions (5.7) and (5.15) can be rewritten as 



E /n^^^'^^^'^^^^P E (w-(7r|g(^*^=(^'»)-W^(7r)' 

i 

Here the statistical weight W{'K\q^'^^^ (r<))) ig defined respecting an additional boundary- 
condition qf(^<^ (^i)) which is imposed in \Ji5^{Ti). We remind that only q^^"^ affects W{'k) 
if TT n 7^ 0. 

The polymer sum in (5.39) describes the interaction between contours F^. It is im- 
portant that every contributing to (5.39) polymer is sufficiently long and contains at least 
[7^3/3] links. In view of (5.30) such a polymer has very small statistical weight. Therefore 
in complete similarity with (5.12)-(5.14) expanding 



and integrating over Yli dq^^ one can derive for (5.39) a representation analogous to 
(5.18). In this representation newly defined diagrams are constructed from contours F 
connected via polymers tt. Taking logarithm and applying Theorem 7.1 one obtains (5.35) 
with correspondingly defined polymers ^. The key fact ensuring the condition (7.3) of 
Theorem 7.1 is that for 7 small enough the sum of absolute values of statistical weights of 
all polymers tt containing given particle and being longer than £3/3 is much smaller than 
the quantity ci2^^'^C entering the Peierls estimate (3.43). 

The route to (5.35) looks rather involved and tedious but it is a standard one in the 
cluster expansion technique. On the other hand we need only minor knowledge about ^, 
namely (5.37). This estimate is obtained by the methods of [DS] and [BKL] along the 
following way. 

The sum over all polymers ^ containing given particle q is equal to 

E E 1^(01+ E E 1^(01 (5-41) 

r:S={r)Bq i-.^BT tt: 7r3g,L(7r)>C(5.3i) ^: SStt 

By (7.6) the first internal sum does not exceed cW'^(r]^\q^^ and the second internal 
sum does not exceed c|VF(7r)|. In turn 

J2 cW''{r]'^\q^^^^^^^)+ c\W{'K)\<-f^'' (5.42) 

r: <5=(r)3g tt: 7r9g,L(7r)>C(5.3i) 

because of (3.37) and (5.31). 

The final result of this subsection can be states now as 
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Lemma 5.2 

The expression (5.3) is equal to 



£'2'-' 
nj 



(5.43) 



i=l 



Expression (5.43) without the last sum is nothing but PJ^'y,p,x{g^^^^q^^^^) defined in (4.21) 
with its main part /3J^'y,/3,a(p^^'^ |g^^'^) given by (4.26)-(4.27'). 

Dobrushin Uniqueness (Basic Calculation). 

Lemma 5.3 

The effective Hamiltonian (5.43) satisfies the Dobrushin uniqueness condition. 

Proof. For the sake of simphcity we first check that the Dobrushin uniqueness condition 
is true for the Hamiltonian J-''y^f3,\{g^^'^^). Then in the next subsection we show that only 
minor modifications are necessary to treat (5.43). 

Take any site x e Xf^ and consider two boundary conditions g^^^ and g^^^ on Z'^^ \ x 
such that they both belong to interval (5.2) and differ only at a site y e Zf^. For the defi- 

niteness we assume that gy^^ > gy^\ Denote by h'^'^^dg^) and i^^^^dgx) conditional Gibbs 
distributions defined by the Hamiltonians J^'yjj,x{Qx\Q^^'^) and ^-■y,p,xiQx\Q^^'^) respectively. 
The Vasserstein distance hetween I'^^^dgx) and I'^^^dgx) is 



dz 



I 



—00 

00 



dz 



r U'\dgx) - iy^'\dgx)) 
J —00 

(i^^'HdQx) - i^^'Kdgx)) 



(5.44) 



where the last equality utilizes the positivity of gx- 

The Dobrushin uniqueness condition is satisfied if one is able to find a function Vxy 
such that 

i?(z.W,z.W)<r,,|4°)-4i)| (5.45) 
and 

(5.46) 



E 



1"xy ^ 1 



To check this condition we follow the strategy of [COPP] and define vtidgx), t € [0, 1] 
as a Gibbs measure corresponding to the interpolated Hamiltonian 



This measure has the density 

exp i - l3J^^^i3^x{Qx\t)) 
P{Qx\t) = — J r , 

Jo exp - (3T^^i3^x{Qx\t)jdgx 



(5.47) 



(5.48) 
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which is difFerentiable in t and therefore from (5.44) we have 

/•I POO pz o 

i?(i.(0),i.(l))< dt dz -^^p{Q,\t)dQ, 

Denoting 



(5.49) 



(5.50) 



one has 

d_ 
dt 



where {•)t denotes the expectation with respect to p(px|0 dQx- Observe that 



(5.52) 



where I^\x, y\Q) is defined by (4.41) (with ^(A \x\y) = ^>(°)(A \x\y) ^ g^'^\A \x\y) 
instead of g) and satisfies (4.42). By direct calculation 



dz 



Qx - {Qx)?jp{Qx\t) dga 



= {Ql)t - (Qx) 



(5.53) 



For sufficiently small 7 and therefore sufficiently large £2 the last expression can be 
estimated by Laplace method and it is equal to 



^2^(0 + 0(^2"''), 



(5.54) 



where ^2^'*(0 ^ value of P-^J-'^^p^x{gx\t) at the point, ^*(t), of minima of 

^'y,l3,x{Qx\t)- As we know from Lemma 4.8 such a minima exists, is unique and lies strictly 
inside interval (5.2). Hence 0{i.2^'^) in (5.54) is uniform in t (see [F]). 
Combining (5.49) - (5.54) we obtain that 

i?(^.(°),^.(i)) < dt{(3g. + 0(^2-'^)) (/f (a;,y) - ^4')(x, y|^)) (^W - ^W) (5.55) 



Setting 



/f (a;,2/) - ^/f (a:,2/|^)|^'dt(/3p.(t) + 0(£2-')) (5.56) 

one concludes that (5.46) is true for sufficiently small 7 because of (3.16), (4.42) and the 
fact that Yliyl'?\xiy) = 1- 
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Dobrushin Uniqueness (General Case). 

We turn now to a complete effective Hamiltonian (5.43) and show that only nonessen- 
tial modification of the calculation of the previous subsection are required to cover the 
complete (5.43). 

The first correction is due to the difference between J-''y,j3,x{Q^^'^^ l^*-^ ■*) and T.y,f3^\{g^^'^^ \ 
q^^"^). This difference is discussed in detail in Lemma 4.6. The corresponding modifications 
of to our previous arguments are the following. 

The difference between volumes of -B^(-) and [B^(-)]^^'^^ produces a factor (1 + 0(7°^)), 

a > 0, in front of •) and ■, ■, •) which clearly is not essential. 

The self interaction like l!^\gx, Qx) or I^\qxi Qxi Qy, Qz) produces nonlinear terms in 
^'y,(5,\{Qx\t) which is defined similarly to (5.47). Their contribution to J^'y,p,\{Qx\t) has the 
form 

C2qI + C^qI + CaqI + t{c5Ql + CQQl){Q^y^ - Q^y^) , (5.57) 

where q does not depend on t, g^; and gy\ Only last, t-dependent, part from (5.57) 
contributes to A^^f^p^x{gx\g^'^\ g^^^) which is defined similarly to (5.50). The corresponding 
effect on dz \Jq -§ip{gx\'t)dgx\ results in terms of the form 



/3c5(4'^ - 4°^) I I (^^ ~ {Ql)t)p{Qx\t) dg. 



and 



13' 



m 



'C6(4'^ - 4°^) dz (^gl - {gl)t)p{gx\t) dg. 



= Pc,{-g^^^ - Q^^^){{Ql)t- {gl)t{g.)t) 

(5.58) 



= /3c6(4'^-4°))((4)t-(^^)t(^.)t) 



addition to £ip(^4'\x,y) - i,4'^(x, y|^)) (^^'^ - g'^^KigDt - (^.)|). Clearly 

\ci\<eir 



(5.59) 



(5.60) 



with some a > 0. 

Applying Laplace method one has 



Qx'^^)t - {Qx)t{Qx)t = 



kg.ity 



+ 0(^2-2'^) 



(5.61) 



where again g^{t) is the minima of J^7,/3,a(^'x|^)- Since 



/3^^7,/3,A(Px|t) =^2P^'+/92c2+/?6c3^a.+/312c4^^+/3t(2c5+6c6Px)(4'^-pf ) (5.62) 

and because of (5.60) all corrections above does not destroy the arguments of the previous 
subsection. 
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As the next step we incorporate in our calculation AHj:^^'^{g\q^^'''i) . According to 
(5.34) this is the same type of a polynomial correction which was just discussed. Thus 
the same arguments work. The necessary smallness, like in (5.60), of coefficients W{D) 
is a consequence of the smallness of the statistical weights of links (5.22)- (5.23) and the 
definition of W{D) (see also (5.30)). 

Finally, to treat X]i=2 3 ^^j^iolQ^^"^) simply observe that for any boundary con- 
dition Q given on \ a; 

i?(iy(^^|^),i/2,3(fo|^)) <7^ (5.63) 

where I'iQxlQ) is the conditional Gibbs distribution given by the whole Hamiltonian (5.43) 
while 1^2,3 (^'x I ^) is a similar distribution given by the Hamiltonian (5.43) without 

5^i=2 3 ^H^^iQl^^^"^) ■ This estimate is obvious in view of definition (5.44) and bounds 
(5.31) and (5.37). One can comment that the contribution of contours or long polymers 
to the free energy of the aiixiliary model is too small to affect anything at all. ■ 



6. Properties of Auxiliary Model 

In this section we use the Dobrushin uniqueness result established for the auxiliary 
model in the previous section to prove Statements 3.3 and 3.5. We begin with a construc- 
tion which is necessary for the proof of both statements. Namely, given a phase a and 
a particle configuration q from the ground state ensemble of the phase a we derive an 
appropriate representation for the logarithm of partition function (3.36). Lemma 5.2 gives 

(6.1) 

Below we consider density configurations, lattice volumes, etc related to the scale £2 and 
in most cases we omit the superscript {£2) from the notations. On few occasions when we 
need other scales we specify them explicitly. 

In Section 4 we studied in detail the minimizers of JP-^,;3,A(^(A)|g(^')) for bounded A 
and A = Zf^. Now we need to study minimizers for 

^7,A(^(A)|g(^^))=/3^,,^,A(^(A)|Q(^^)) + Ai/«(KA)|Q(^^)) (6.2) 

where both and ^(A) belong to the ground state ensemble of the phase a. This 

is a finite range translation- invariant Hamiltonian similar to J^^^p^x{g{A)\q^^ ^) but not 
allowing a simple representation of the type (4.35). For a bounded A the existence of 
at least one minimizer of Hj^\{g{A)\q^^ ^) follows from the compactness of its domain, 
{Q/3,(t — C7 Qf3,(T + C)^- This domain is convex and the minimizer is unique because the 
function Hj^x{g{A)\q^^''')) (of the finite number of variables g^-, x e A) is also convex. This 
is a consequence of estimate (3.16) from which it is not hard to see that for any q^^"") and 
g{A) the Hessian matrix of H^^\{-\q^^''^) calculated at g{A) is positive definite with the 
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mass bouded from independently of ^(^"^ and q(A). Depending on the context we denote 
the unique minimizer of H^^x(g(A)\q^^'''>) by g, q{A) or 

To clarify the structure and convexity of H-y^\{g{A)\q^-^''^) we introduce 

AH,MA)\q^''^^) = H,MA)\q^''^^) - H^,xmW^^) (6-3) 
It is not hard to check that 



oo ^ 



+4 Yi ^(^) E n ^-(^) n 



(6.4) 



Here A^^; = gx — Qx{A) and D is a connected family of links with endpoints forming a set 
X[D) C Ui^. The number of links in D is less than C(5.3i). The subsets Xi[D) and X2{D) 
form a partition of X{D) with at least two elements in X2{D). The second internal sum 
is taken over all partitions of that type. For W{D) one has the estimate 

E \WmQ^^^^^^<l'^ (6.5) 

D: X(D)3x 

as follows from (5.30). Also 

|A^,|<2C (6.6) 

The quadratic part of (6.4) is 

^4 J2 ^^^'-^^2 E 7(^n^i,^2|^(A))A^.,A^., (6.7) 

with 

,-4tt -E^^'H^i,^2|^(A)) > mi/3) > (6.8) 

for any xi e [A]*^^^)^ jjere /(^^(a;i, a;2|^(A)) is an analogue of (4.41) and satisfies the same es- 
timate (4.42) for 7 small enough. ( By construction |/(^)(a;i, X2|^(A)) — /*^^)(xi, X2|^(A))| < 
c7"2/(2)(2;^^ 2:2) ). Moreover, in view of (6.5) and assuming that ^ is chosen to be small 
enough with respect to m(/3) 

A#,,,(^(A)|g(^^)) = Y ^^^l + U{g{A),^Qm (6-9) 

a;e[Al(^2) 

with positive convex U{g{A), Ag{A)) having minima at A.g{A) = 0. 
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Along with the minimizers of H^^x{g{A)\q^^''^) for bounded regions A we also need 
a global minimizer, i.e. an analogue of Qis^cr- To find such a minimizer and to show its 
uniqueness ( for given a ) we consider a suffuciently large (with respect to the range of 
Hj,x{-)) periodic box A and the corresponding Hamilltonian Hj^x{q{A.)). All convexity 
considerations remain true for H^^x{q{-^)) such that Hj^xiQi-^)) ^ unique minimizer 
^(A). Because of the translation invariance of H^^xiQi-^)) this minimizer is a constant 
configuration Qxi-A) = g, x & A. (Otherwise one is able to construt another minimizers by 
space translations). Since -ff-y a(^(A)) is of finite range we observe that g is independent 
on A for all sufficiently large periodic A and hence g is the global minimizer we are looking 
for. Note that the specific energy r^x{s) = lim^^^^^d \A\~'^ x{q{^)) of any constant 

density configuration ^(A) = s is just a finite sum and g is nothing but the value of 
s at which h^^p^xis) achieves its minima. Thus given a such a minima is unique for 
s e iQf3,(j ~ CiQp,(T + C) cind we have different values g/3,x,- and g)3,x,+ corresponding to 
cr = — 1 and a — +1. 

Denote by A(/9, 7) the value of A at which hj^i3^x{Qp,x,-) — h^,f3,\{Q(3,\,+) and set 
Ql3,a — Qp x{f3) a- "^^^ important consequences of the representation of g/3,x,a via h^,p,x{s) 
are the existence of A(/3, 7) and the estimates 

5.7.1 \gi3,a - Qp,X,a\ < C7"^/' (6.10) 

and 

5.7.2 |A(/3)-A(/3,7)|<c7"^ (6.11) 

Indeed, let h^^j3^x{s) be the specific energy oi g = s calculated via the Hamiltonian 
^7,a(^(A)) + /\H^^l{g{A)) with periodic A. Then 

5.7.3 \h^,i3,x{^)-h^,p,x{s)\<ci^'' (6.12) 

as follows from (5.31). 

Another way to calculate h~^,j3,x{s) is the approach of Section 4. Here the starting 
point is the partition function (5.15) with fixed density configuration g'^^'^^ = s. Then 
one aproximate this partition function by using the density configurations defined with 
respect to the scale £\ <^ £2- (Note that the density configuration which is constant in 
the scale £2 is not necessarily a constant one in the finer scale £1). As the first step of the 
approximation one shifts all particles to the centers of the corresponding boxes C*-^^-' for 
the price of the error (in the value of H^,x{q^^^\^)) + /^H''^^{g^^^\A))) which in absolute 
value does not exceed C7"^|A| (see Lemma 4.5). Lemmas 4.6 and 4.7 then show that up 
to the same error the Hamiltonian H^^x{q'^^^\A)) + /^H^^'^^{g'^^^\A)) can be approximated 
by (4.26). 

Now the convexity considerations (see (3.5)) imply that among all density configura- 
tions ^'^^^-'(A) such that ^'^^^-'(^'-^^-'(A)) = s the configuration g^^^\A) = s has the minimal 
value of the mean field functional (4.26). Hence (see (4.19) and above) counting only the 
contribution of g^^'^^A) = s produces the error which again in absolute value does not 
exceed 07^*1 |A| (in fact it is much smaller). The specific energy of the constant density 
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configuration g^^^^ = s calculated via (4.26) is nothing but the mean field specific energy 
-P/3,a(s) and therefore 

\F^A^)-h^,p,x{s)\<cr' (6.13) 

Since /i7,^,a(s) = Fp^x(s)+^^^2 ^kS^ with |a/s| < 07°^^, i.e. h^,p,\{s) as a small polynomial 
perturbation of Fp^s), and 

\F^,x{s)-h^,0,x{s)\<cr\ (6.14) 

as follows from (6.12) and (6.13), one immediately obtains (6.10) and (6.11). ( The exis- 
tence of A(/3,7) follows from an elementary linear analysis and continuity of /i-y,/3,A('S) )■ 

Finally observe that a straightforward modification of the arguments of Lemma 4.8 
shows that \gi3,x,a — q'^ \ satisfy (4.36). 

With gf3,\,a and g'^ being properly defined we construct now a convenient repre- 
sentation for the log of the partition function (6.1) using the following interpolation trick. 
For t e [0, 1] and s e [0, 1] introduce an interpolated Hamiltonian 



00 ^ 

Ai?,,,(A,(A)|,-(^=);M) = ^^ E n,---^ ^"- 



+tei Yi ^(^) E n ^-(^) n (6.15) 

+ .E^<A(KA)|g(^^)) 



i=2,3 



The Hamiltonian A.H^AAg{K)\q^^ ^]t-,s) satisfies the Dobrushin uniqueness condition 
similarly to Ai7^,A(A^>(A)|^(^'); 1, 1). We denote by (■) 

t s A qC^") expectation with 
respect to the corresponding Gibbs measure in the finite domain A with the boundary 
condition q^^"^ and we denote by {■)t,s the expectation with respect to the corresponding 
unique limit Gibbs distribution. 

From (6.1) and (6.15) by direct calculation one obtains 

logZ^^^,,(A|g(^^)) = -^^,,(^(A)|g(^^)) 

+ 4 E ^(^) E n n ^^-)t,o,A,.-(-) 

D!Z:[A=](^2) x€Xi{D) ''^ x€X2{D) 

(6.16) 



+ E W{t) f\s{P{T,Ag)),^, 
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Here ^(A) = is the minimizer of iy^,A(-|9^^"^) and notations P(r, A^) and W{^, Ag) 

are used instead of -P(t) and W{^) to underline the density configuration ^(A) + A^(A) 
with respect to which these quantities are calculated. We use two parameters s and t 
instead of a single one for the technical transparency. In models {■)t,o,A,q('^''h i-e- with 
s = 0, the range of interaction is finite so it is simpler to control decay of correlations. 

Namely, a standard consequence of the Dobrushin uniqueness theorem [D2] is the 
estimate 

CO n—1 

{Ag,)t,o\<2CYl H n^f^f^+i' (6.17) 

n=2 i=l 

H2--.!/„-ie[A](^2) 
yne[Ac](^2) 



where rxy are defined by (5.56) with nonessential modification discussed below (5.62). 
Similarly 



n-l 



< (2C)|"| E E n 

n=2 vi&x i=i 

M-.-.J/n-ielAjC^a) 

Hne[AC](^2) 



ViVi + l 



(6.18) 



Now it follows from definition (5.56) of Vxy and estimate (3.16) that 



, / A \ / A \ , cC / 1 + a{j3) \ 

l(A^a.)t,0,A,,-(A=) - (A^a.)t,0| < 2 ) 



[7dist (a;,A=)] 



(6.19) 



and 



[Tdist (X,A=)] 



(6.20) 



which are similar to (4.36). Also given x 



n-l 



E E n 

n=2 !'i=^ i=l 

J/2 ■•■.!/« 62^2 



l-a{p) 



(6.21) 



In fact it easily follows from the structure and smallness of |VF(t)| and |VF(OI 
pressed by (5.31) and (5.37) that the same exponential decay of correlations takes place 
for models with s ^ which include the infinite range part s X]i=2 3 
the Hamiltonian (6.15) 
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The final representation of logZ;^^ ^(A|q'(^'')) is 



logZ^\,(A|^(^^)) = -i?^,,(^(A)|g-(^^)) 

-ii Yi ^(^) E n n 

- E W{t) l'\s{P{T,Ag)),,, 

T0A=: -L(r)>c ° 



x€X2{D) 



T0A<=: L(r)>c ° 

5.14 (6.22) 
Here speaking about (/(A^ia;))^ ^^^^^(ac) we assume that Ag^ = gx — fe(A) while for 
{f{Agx))t,s we mean A^i^: = Qx - Qi3,\,a such that {f{Agx))t,s is truly independent on 
any boundary conditions. 



Proof of Statement 3.3. From the interpolation trick we have a representation for the 
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metastable free energy 



fa,\n = ^2^/3,A,-T(log^/3,A,a - 1) - ^2A^/3,A,a 



■2\ 



'^2 



D^I,(2),I,(4) 

5.15 r:^r)^o\^^^>\ (6.23) 



|X2(I3)|>2 X2(I3) 



r: X(t)30 

y: /'d. ^'^'^ 

?:X(030 







1^(01 



Here the definition of X{t) and X{l^) is similar to that of X{D) and is defined with 
respect to Qi3,\,a- 

The difference ^^(A|g(^')) between RHS of (6.22) and /^a.iI^I 

can be estimated as 
|^a(A|g(^^))| < c^^^^£^7-"^ < c7->n A^l (6.24) 

•^2 

Here |AnA'^| is the (d— l)-dimcnsional volume (area) of hypersurface separating A and A'^ 
and |AnA''|£2 '^'^S""^ is the number of lattice points in [5N\^^'^\ The first factor cl^ is the 
upper estimate for the sum of the absolute values of the error terms associated with given 
site X e [dA]^^^\ To be more precise, the contribution to ga{A.\q^^''^) is given by terms in 
(6.22) and (6.23) crossing the boundary of A directly or indirectly. Directly crossing terms 
are D-s with X{D) intersecting both A and A^, r-s with X{t) intersecting both A and 
A"^, ^-s with X{^) intersecting both A and A"^ and links with endpoints in both A and A^. 
Indirectly intersection terms come from estimate (6.18) for the difference terms in the last 
three lines of (6.22). Graphically those terms can be represented by chains of r-s joining 
something inside A, say some D e A, with the sites outside A. It is important for us that 



52 



geometrically any of the terms crossing the boundary passes through some site x e dA. 
This site x may be: one of the endpoints of the link, contributing itself or as a part of 
r or ^, or the endpoint of some r^y or it may belong to the support of some contour 
which is a part of some polymer ^. For all these involved weighted objects we always keep 
the property that the sum of the absolute values of the statistical weights of all objects 
passing through given point is less than absolute constant. Multiplying this constant by 
the factor which depending on the notations is present explicitly or implicitly in front 
of the sums just discussed, we reproduce (6.24). We need now some additional work to 
improve on (6.24). 

The simplest observation is that using (5.43) and (5.37) one obtains that the sum over 
the objects containing r with L{t) > 0(5.31) or ^ as a constituting element is less than 07^^^. 

Another source of smallness are various expectations {Y[^ ^Qx) in the corresponding 
terms. In particular we check that given x and for any k > 1 

{\Ag,\'^) < c(£^)-'=A (6.25) 

where the expectation means any one of those contributing to (6.22) and (6.23). The 
power — ^ is taken for the definiteness only and it can be replaced by — | + e. We prove 
(6.25) adapting to our situation relatively abstract Lemma 8.1 in the spirit of [R2]. The 
correspondence between current notations and those of Section 8 is given by the following 
list of analogous objects: (f)^ and A^^;, k and 0. and 2^, M and 2(£2j R and C(5 31)7"^, J^y 

and gl^'^^x^ylo) and so on. The bound (6.25) is an easy consequence of Lemma 8.1 (which 
deals the finite range interactions) as the infinite radius part, s ^^=2 3^Hj\iQ^^'^''\'l^^ '*)> 
of (6.15) is so small that it can not increase the expectation of |A^j;| more than in 1 + 07^'' 
times. Note also that in (6.22) and (6.23) we have terms (|A^a;|'^) only with k >2. 

In a similar way applying estimate (8.12) to \ {lla:eX2{D) ^Qx)t,o,A,q^''"^~ 
x.,{D) ^Px)t,o| we obtain improved (6.20). Indeed, the RHS of (6.18) is nothing but 
the upper bound for the Vasserstein distance and a square root of this estimate produces 
an exponential bound just with twice smaller exponent than in (6.20). 

These improvements reduce by the same factor, (^2)"^^) ^^e estimates of the cross 
boundary terms coming from lines 2 and 5 in (6.22) and line 7 in (6.23). 

As an immediate consequence of these improvements we can see that that the con- 
tribution, e(A), to the difference £2'^if+ x 7 ~ /- a 7) coming from the non energy terms ( 
lines 7,8,9 in (6.23) ) is a continuous function of A and |e(A)| < 0(^2)"^^ • The energy parts 
( lines 1-6 in (6.23) ) of ^ and £2'^/^^ 7 coincide with each other at A = A(/9,7) 

and they contain explicitly terms which are linear in A. Hence shiftig A(/3, 7) by at most 
0(^2)"^^ one can find the solution, A = \{(3,'y) of (3.41) which proves Statement 3.3. ■ 

A useful consequence of our improvements is the estimate 

^,(A|g(^^)) + (^^,a(^/3,A,7(A)|^/3,A,7(A'=)) - ^[/7(^^,A,7(A)|^/3,A,7(A^)) 

-^^,A(^(A)|g(^^^)) <c7-^(£^)-2A|AnA^| (6-26) 
<7i/2|AnA'=| 
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where Uj{Qi3^x^j{A)\gf3^x^j{A^)) is defined as in (2.11) but in terms of the Hamiltonian 
H{-). Observe that Hj,\{g3,\^-y{A)\gf3^x^y{A^)) — H^^x{g{A)\q^^'''>) is rather small while the 
compensating term —^U^{Qi3^x,-y{A)\Qi3,x,'y{A'')) is of order C7~-'^|A n A*^]. 

Proof of Statement 3.5. To check estimate (3.44) consider Int^(r) = 
Um: <T^(r)/a(r)Intm(r) and introduce a strip S{T) = d^^^/^hnt^ {T) \ S^{T), where 6^{T) 
is defined similarly to (^^(F) (see (4.64) and below). First we rewrite the numerator of 
(3.38) as 



dqTL , e-/'^^.^(^l^-^^^^^»)Z,^,,,(lnt^(r) U5(r)|g(^"(r))) (6.27) 

Q(Supp(r)\s(r)) V{<1)=V 



Observe that for any q contributing to (6.27) the restriction of q to S^{T) belongs to the 
ground state ensemble of the phase —a{T). Given such q(^^(^)) denote by ^(Int^(r) U 

S{T)) = gi'" the minimizer of (6.2). Then using representation (6.22) and definition 
(6.23) one can see that 

Z^,,A^nt^i^)^S{T)\q^'"^^^^)^Z^^^^^^^ 

xexpfc'^/4(r),A,-y|Int^(r)| 



-it/(^^^^"-^(int^(r))|^^'^"-'(aint^(r))) ^^^^^ 

-(#„.(^^'^"'^»(int^(r))|^^'^"'^»(aint^(r))) 
-iy7,A(^^,A,.(int^(r))|^^,A,.(5int^(r)))) 

+ ^_.(r)(Int^(r)|^^^'^''»(aint^(r)))) 

The term ^_o-(r)(Int^(r)|^^^'' ^^'^ (5Int^(r))) is the sum over "polymer type" terms cross 
the boundary of Int^(r) and it can be estimated (see (6.26)) by 

|^-a(r)(Int^(r)|^^'*^''^"(aint^(r)))| < C7i/2|lnt^(r) n Int^(r)'^| (6.29) 



On the other hand the strip S{T) is so wide that in the 7 ^ neighborhood of Int^(r) fl 

Int^(r)'= the difference |^^''' '"^^^ — gf3,x,a\ is exponentially small as follows from a straight- 
forward analogue of (4.36). Thus 

</3,A(lnt^(r) U 5(r)|g(^^(^))) = Z,V(^(nk^'^^^^\ ^/3A-.(In^^ 

xexp(^2"'/-.(r),A,7|I^t''(r)l 

1 (6.30) 

--Uigp,xA^nt^(mQp,xAdlnt^im 

+ ^_.(r)(Int^(r))) 
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where 

5.36 |^_^(r)(Int^(r))| < c7^/2|int7^(P)Ulnt^(r)| (6.31) 

Here ^_^(r)(Int^(r)) collects the contribution of ^_^(r) (Int^(r) |^«^'''''^^' (aint^(r))) to- 
gether with corrections due to replacement of U{g'^^^ "^"^'^ (Int^(r)) |^^''' "^"^^ (9Int^(r))) by 
t/(p/3,A,a(Int^(r))|^;3,A,.(aint^(r))) and the estimate of (Int^(r))| 

^'^''^''»(aint^(r))) - ^^,,(^^,,,.(int^(r))|^;3,A,a(aint^(r))) 

Using (6.10) we obtain from (6.30) our last estimate 

xexpfc'^/4(r),,^^|Int^(r)| 

5.37 ^ (6.32) 

--c/(p^,,,.(int^(r))|^^,,,.(aint^(r))) 

+ ^-a(r)(Int^(r)) + A 

where |A| < c^^/^\S(r)\. 

A similar estimate is true for the denominator of (3.38) which after integration over 
q{S^{^) implies (3.44). ■ 



7. Appendix. Polymer Expansion Theorem 

Consider a finite or countable set the elements of which are called (abstract) dia- 
grams and denoted 9, 9', ets. Fix some refiexive and symmetric relation on x 0. A pair 
6*, 6*' G X is called incompatible {9 7^ 9') if it satisfies given relation and compatible 
{9 ~ 9') in the opposite case. A collection {9j} is called a compatible collection of diagrams 
if any two its elements are compatible. Every diagram 9 is assigned a complex-valued sta- 
tistical weight denoted by w{9), and for any finite A C an (abstract) partition function 
is defined as 

A.1 z{A)= Yi n«'(^^)' (7-1) 

{e,}cA J 

where the sum is extended to all compatible collections of diagrams 9i G A. The empty 
collection is compatible by definition, and it is included in Z{A) with statistical weight 1. 

A polymer tt = [6*^'] is an (unordered) finite collection of different diagrams 9i & Q 
taken with positive integer multiplicities e^, such that for every pair 9', 9" G tt there exists 
a sequence 9' = 9i^, 9i^,. . . , 9i^ = 9" E tt with 9i. 7^ ^ij+i > J = 1) 2, . . . , s — 1. The notation 
TT C A means that 9i E A for every G tt. 

With every polymer n we associate an (abstract) graph G{n) which consists of J2i 
vertices labeled by the diagrams from tt and edges joining every two vertices labeled by 
incompatible diagrams. It follows from the definition of G{7r) that it is connected and we 
denote by r(7r) the quantity 

r(7r) = n(^,!)-i E (-1)""'' (7.2) 

i G'CG{tv) 
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where the sum is taken over all connected subgraphs G' of G(7r) containing all of X^jSi 
vertices and denotes the number of edges in G' . For any ^ e tt we denote by e{6,7T) 
the multiplicity of 9 in the polymer tt. 

The polymer expansion theorem below is a modification of results of [Se] and [KP] 
proven in [MSu]. See also [D3] for similar results. 

Theorem 7.1 Suppose that there exists a function a{9) : 6 i— > R""" such that for 
any diagram $ 

\w{e')\e''^^'^ <a{e). (7.3) 

Then, for any finite A, 

logZ(A)= 5^«;(7r), (7.4) 

ttCA 

where the statistical weight of a polymer tt = [9^^] equals 

w{tt) = r(Tr)Y[w(9iy\ (7.5) 

i 

Moreover, the series (7.4) for logZ'(A) is absolutely convergent in view of the estimate 

J2 e{9,n)\w{7T)\<\w{9)\e'^^'\ (7.6) 

tt: ttBO 

which is true for any diagram 9. 

Corollary 7.2 For any function b{9) : E+ consider modified statistical weights 
of diagrams 

w{9) ^ w{9)e^^^\ (7.7) 

and suppose that still 

J2 \w{9')\e''^^'^ <a{9). (7.8) 

6': d'i^e 

Then for any family, II, of polymers such that any tt e 11 contains given diagram 9 



e{9,7r)\w{7r)\ < \w{9)\e''^^^ mm I FT e^^^'^ | (7.9) 

7r€n \ / 



Tren \6>'e7r 
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8. Appendix. A Technical Lemma 

Consider a spin model on the lattice given by the formal Hamiltonian 

m) = E '^^ + ^ E ^-y'^-'^y + E ^^'^^ (8.1) 

Here the spin variable takes M discrete values, including 0, from the bounded interval 
[—a, a]. The last sum runs over some family A of sets A = {xi, . . . , X\j^\ G Z^} containing 
l^l not necessarily different sites x and (f)A — YIxea'^x- The interaction, J^y and Ka, is 
of finite range i? < oo, i.e. Jxy = if dist(a;, y) > R and Ka = if diam(^) > R. Suppose 
that 

^\Ja;y\ = l-a, 0<q;<1 (8.2) 

and 

E^'^'l^AKa (8.3) 

A3x 

and denote by //(0(A)) the corresponding Gibbs distribution in the finite domain A C Z'^ 
with zero (= empty) boundary condition in A*^. 

Lemma 8.1 

For any c > 1 and sufEciently large inverse temperature 

M + R 

K > (8.4) 

c 

one has 

At(|0x| > «"^/'^) < 6"'='^'^', \/xeA (8.5) 



Proof. For an arbitrary configuration (f>{A) define spots S{(j){A)) as i?-connected compo- 
nents of sites X E A with |0a;| > Taking x E S set 

h{x, 0(A)) = E Jxy(l>x(t>y'S... , , + E |>„.n U I (8-^) 
y^x A^A: ABx 

Then 

i? (0(5)|0(A \ <5)) = E m) (8.7) 

x£S 



and 



h{x, 0(A)) > (8 8) 



Here in (8.7)- (8.8) we used (8.2), (8.3) and the fact that |0a;| > \(f>y\ for any x E S and 
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Now the probability that 0(A) contains a spot S with the fixed value of (f){S) can be 
estimated as 

,,,,,, ^HA\s) exp ( - KH{cf>{SmA \ S)) - kH{^{K \ S))) 

E<^(A)exp(-«iy(0(A))) 

< cxp ( y """^ ^ ~ ^"^"^^^ ^ ^^^^ 

5.26 2 ^ " E0(A): <^(S)^0 exp ( - ((/-(A))) (g.Q) 

a 1/6, „,x E<^(A\5)exp(-Ki?(</'(A\^))) 



< exp ( - ^«V6|5|) 
= exp(-|KV6|5|) 



E0(A): 0(S)=O exp ( - KH{(j){A))) 



Hence the probability of the spot S 

5.27 i,{S) = l^i^iS)) < (c(8.4)^)l^' exp ( - l^i/^l^l) < exp ( - ^k'/'\S\), (8.10) 

where the last inequality is true for k large enough. 
Finally for any x E A 

MI</>.|>^-'/^') = 5^M^) 

SBx 

oo 

< E(^(8-4)-)'"exp(-^«^/^|^|) 

5.28 |5|=1 , (8.11) 

oo 

< ^exp(-V/6|^|) 

|S|=1 

< e '^'^ 

where again k is large enough. ■ 

Let now ijS^\4>xi 4>y) ^md ijS^\4>xi 4>y) be a distribution at sites x, y e A of conditional 
Gibbs measures in A with two arbitrary boundary conditions ^'^'^^(A'^) and (()'^^\h.'^) re- 
spectively. The Vasserstein distance R{pi,'^^\ ^'^^^) between ^'^^^ and /U*^^-* is given by some 

coupling fi{(f)^\ (^y^-*; (f)^\ (j)^^)- Denote by {■)^^\ i')'"'^'^ cind (•) the expectation with respect 
to fi^^^ and n and suppose that 2a < 1. Then 

. x\ 1/2 / . .\ 1/2^ 

xmax (i?(/.W(0,),/.W(0,))V2,i?(^(o)(<^^),^(i)(<^^))i/2), 

5.29 (8.12) 
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< max 



Similarly one can estimate via the Vasserstein distance and the second moment the differ- 
ence KH^^x <^.)^°^ - iU.ex any X C A. 
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